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Using a recently suggested method of bosonization in an arbitrary dimension, we study the anoma- 
lous contribution of the low energy spin and charge excitations to thermodynamic quantities of a 
two-dimensional (2D) Fermi Uquid. The method is slightly modified for the present purpose such 
that the effective supersymmetric action no longer contains the high energy degrees of freedom but 
still accounts for effects of the finite curvature of the Fermi surface. Calculating the anomalous 
contribution Sc(T) to the specific heat, we show that the leading logarithmic in temperature cor- 
rections to 5c{T) /T^ can be obtained in a scheme combining a summation of ladder diagrams and 
renormalization group equations. The final result is represented as the sum of two separate terms 
that can be interpreted as coming from singlet and triplet superconducting excitations. The latter 
may diverge in certain regions of the couphng constants, which should correspond to the formation 
of triplet Cooper pairs. 

PACS numbers: 71.10.Ca, 71.10.Ay, 71.10.Pm 



I. INTRODUCTION 

At low temperatures, thermodynamic properties of 
fermions with a repulsive interaction bear strong re- 
semblance to those of an ideal Fermi gas. This is the 
quintessence of the Landau theory of the Fermi liquid^. 
It is assumed in this theory that interaction effects merely 
renormahze quantities such as the fermion mass or the 
density of states. In fact, this renormahzation can be 
large for a strong interaction, thus making perturbative 
methods inapphcable. However, such obstacles are al- 
ways overcome once the renormahzed quantities are re- 
placed by phcnomcnological effective parameters. 

Fohowing the similarity to the ideal Fermi gas, one 
could expect that quantities like 

5(T) = ^, (1.1) 

where c(T) is the Fermi liquid's specific heat, or its spin 
susceptibility x{T) had to be analytic functions of r^/e|, 
with ep being the Fermi energy and T the tempera- 
ture. [In the leading order in T, one should have a fi- 
nite g (0) considered as a phenomenological Fermi liquid 
parameter.] However, several studies revealed for three- 
dimensional Fermi liquids the existence of corrections to 
the specific heat c(T) of the order T'^lnT, which is in- 
compatible with ideal Fermi gases i^"— In three dimen- 
sions, there are also logarithmic contributions |qpln|q| 
to the non- homogeneous spin susceptibility xi"^-, where 
q is the wave vectori^ 

In two dimensions (2D), non-analytic corrections to 
the quantities g[T), Eq. (jl.ip . and x{T) are stronger 
and, in the lowest order in interaction, have been found 
to be proportional to T.— These anomalous contribu- 
tions were attribute d^^d'^ to one-dimensional backscat- 
tering processes imbedded in the two-dimensional mo- 



mentum space. The linear in T correction to g{T) has 
been verified experimentally in a '^He fluid monolayer 

The problem of evaluating the anomalous contribu- 
tions was reconsidered in Ref. [l^ with the help of a 
supersymmetric field theory especially designed to de- 
scribe those low energy bosonic excitations which are re- 
sponsible for the anomalous contributions to the ther- 
modynamics. It was found that earlier calculations^r— 
had not been complete and so far unforeseen logarith- 
mic contributions to 5c{T)/T^^ were discovered for di- 
mensions 0? = 2, 3. Similar anomalous contributions have 
been found for the spin susceptibility x(T) using either 
the supersymmetric method mentioned above^^ or the 
conventional diagrammatic technique-^. Both the meth- 
ods led to identical results for the spin susceptibility in 
2D models. 

In dimension d = 1, the supersymmetric approach to 
find the spin susceptibilitj*i^ reproduced the results of 
earlier theoretical work o^^'^° . As to the specific heat of 
a one-dimensional Fermi gas, the result of the supersym- 
metric field theory of Ref. [l5| can be mapped on known 
results for the Kondo modeU^ or for the XX Z spin-i 
chain^Sj showing agreement. A more recent study^ con- 
firmed the supersymmctry approach using the conven- 
tional diagrammatic technique. 

In spite of tlie agreement between the results obtained 
by these different methods in one dimension, a direct di- 
agrammatic computation of the anomalous specific heat 
carried out up to the third order in the fermion-fermion 
interaction by Chubukov and Maslov (CM) for the 2D 
Fermi liquid^^ led to a result that did not coincide witli 
the one obtained by the supersymmetric approach in 
Ref. [l^. Both of them contained logarithmic corrections 
to the anomalous contributions. However, while in the 
framework of the supersymmetry method of Ref. [Tsl . the 
non-trivial contributions originated purely from spin ex- 
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citations, CM obtained contributions from both spin and 
charge excitations. They attributed the difference be- 
tween the resuhs to the fact that not all effects of the 
finite curvature of the Fermi surface had been properly 
taken into account in the approach of Ref. [isl 

Of course, curvature effects are absent in one dimen- 
sion and charge excitations do not influence the spin sus- 
ceptibility xiT) in any dimension in the non- logarithmic 
lowest order in the interaction. As a result, no discrep- 
ancy could be seen in these cases. Nevertheless, since ev- 
idently certain effects of the finite curvature of the Fermi 
surface were neglected in the supersymmetric method of 
Ref. [l^ it is important to find the correct way of cal- 
culations. At the same time, approaches based of the 
conventional diagrammatic expansions for fermions be- 
come inapplicable beyond some low orders of perturba- 
tion theory. Indeed, CM performed a standard perturba- 
tion theory to third order and treated higher orders by 
plausibility arguments 

In this paper, we revise the approach of Ref. [l^ for the 
2D Fermi liquid taking into account all the necessary ef- 
fects of the curvature of the Fermi surface. As a result, we 
are able to sum up all leading logarithms, thus correct- 
ing the previous result for the function g{T), Eq. (|l.ip . 
To third order in the interaction potential, our result 
agrees with the conventional perturbative calculation^^. 
Moreover, our result in all orders in the large loga- 
rithm ln(ej;-/T) shares the same asymptotic behavior as 
the conjecture suggested by CM. In principle, our method 
and results are applicable for both repulsion and attrac- 
tion unless one reaches a singularity in the final formulas. 
We argue that the singularities, if existing, correspond 
to the singlet or triplet Cooper superconducting pairing. 
Remarkably, the final formula for the function g (T) con- 
tains a sum of separated spin singlet and spin triplet 
excitations. It is important to emphasize that the mod- 
ification concerns the dimensions d > 1 only, whereas 
the method leads for d = 1 to the same results as those 
obtained in Ref. [T^. 

The calculations are performed using a modification of 
a recently suggested bosonization scheme of Refs. [23ll23 . 
In contrast to these previous works, we derive an ef- 
fective supersymmetric action describing only low-lying 
modes. This is achieved by singling out the slowly vary- 
ing pairs of the fermionic field in the interaction term. 
Subsequently, we decouple this interaction by means of 
Hubbard- Stratonovich auxiliary fields slowly varying in 
space and imaginary time — similarly to what was done 
in Ref. [l^. Here, however, this decoupling is followed by 
the derivation of equations of motion using the method 
of Refs. [23II23 . In contrast to the equations of Ref. [H, 
the present equations preserve all necessary effects of the 
curvature of the Fermi surface. 

The solution of the equations of motion is represented 
in a form of an integral over superfields "I* , which do not 
only depend on conventional coordinates r and imaginary 
time T but also on anticommuting variables 6,0*. This 
integral representation allows to average over the auxil- 



iary fields before we obtain the final effective field theory 
for the low energy bosonic charge and spin excitations. 
Such a representation, suggested in Refs. [23ll2^ differs 
from the supervector representation used in Ref. [l^ and 
is considerably more convenient for explicit calculations. 

Although the general calculational scheme based on 
this superfield action shares certain similarities with that 
of Ref. [THI, the finite curvature of the Fermi surface sup- 
presses several otherwise logarithmic contributions. Con- 
sequently, different final results are obtained as a result 
of a different calculational procedure. For instance, the 
quartic part of the action can be renormalized by sum- 
ming ladder diagram series instead of solving rcnormal- 
ization group equations. 

The calculations performed here can be important not 
only from the point of view of finding the complete pic- 
ture about the anomalous contributions to the thermo- 
dynamics of the 2D Fermi liquid but also as a demonstra- 
tion of how the higher-dimensional bosonization scheme 
suggested in Refs. [231123 can be used as a method in an- 
alytical studies. The experience gained on this compara- 
bly simple example may become important for attacking 
more difficult and more interesting problems of strongly 
correlated systems. 

The paper is organized as follows: In Sec. [TTl we de- 
rive the effective low energy field theory for the anoma- 
lous thermodynamic contributions. Starting from a gen- 
eral model of repulsive interaction, we discuss and single 
out the relevant soft modes and bosonize the microscopic 
fermion model in the low energy limit. 

Section [Hi] discusses the leading perturbative correc- 
tions to both the thermodynamic potential and the ver- 
tices of the low energy field theory on one-loop level. We 
identify the logarithmically divergent one-loop diagrams 
that are important for the subsequent renormalization 
group analysis. This analysis is presented in Sec. IIVI 
in which we derive and solve the flow equations for the 
coupling constants of the low energy field theory. 

In Sec. |Vl we apply the bosonic technique and the 
results of the renormalization group analysis to evalu- 
ate the anomalous contribution to the specific heat be- 
yond the T^-term obtained from second order perturba- 
tion theory. First performing an explicit perturbation 
expansion to third order in order to check once more our 
bosonic approach, we eventually include the completely 
renormalized vertices and find the non-analytic contribu- 
tion to the specific heat in all orders in \n{eF/T). 

Concluding remarks are found in Sec. I VII 



II. LOW ENERGY FIELD THEORY 

In this section, we formulate the microscopic model 
for the interacting fermions and derive the low energy 
field theory that catches the non-trivial physics of the 
low-lying bosonic excitations. The derivation does not 
require to specify the dimension d of the system and we 
assume d to be arbitrary here. 
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A. Microscopic fermion model 

We consider a gas of spin-i fermions described by the 
Hamiltonian 

H = Ho + Hint, (2.1) 
where Hq is the kinetic energy, 

Ho = J2[ 4(r) [e (-iVr) - /i] c.(r) d^r . (2.2) 

In Eq. (12. 2p . r and cr = ±1 denote the coordinates and 
spin, respectively, (r) [co- (r)] are creation (annihila- 
tion) field operators, and jjL is the chemical potential. 
In the simplest case, e (p) = p^/277i with m being the 
fermion mass. In this case, the Fermi surface is a (d— 1)- 
dimensional sphere. For a more general spectrum e (p) , 
the Fermi surface has a more complex shape but this 
does not lead to a qualitatively different physical picture 
as long as the Fermi surface remains smooth and there is 
no nesting. 

The second term in Eq. (j2.ip stands for the fermion- 
fermion interaction and takes the standard form: 

i?i„t = \Y.j cl{v)cl{v')V{v~v')c,.{v')c,{v) d-^rd^r' . 

(2.3) 

At the moment, we do not specify the form of the function 
(r — r') except for its positivity, guaranteeing repulsive 
interaction. 

Equations (j2.ip - (|2.3p constitute the model in the 
Hamiltonian form. It is more convenient for our purposes 
to use a functional integral representation, in which the 
partition function Z is written as 

Z = j cxp{-5o -5int} V{x\x) ■ (2.4) 
Herein, the Euclidean action is given by 

x[dr+e (-iVr) - fi] Xa{v, t) A'^vAt, (2.5) 

'^int = ^ V / / Xa{Y, t)xI> (r', t)x^' (r', T)Xa{v, t) 

X V{y - r') d'^rd'^r'dr . (2.6) 

In Eqs. (|2.4|) " ()2.6|) . P = 1/T is the inverse temperature 
and XiX* ^-re Grassmann fields which are antiperiodic in 
imaginary time r, x(t + /3) = — x(''')- 

Equations p.4p - p.6|) are the starting point for our 
analysis. 

B. Low-lying modes 

The spin and charge excitations at low temperatures, 
which we are interested in, correspond to the low-lying 




FIG. 1: Decomposition of the interaction <Sint, Eq. (|2.6|l . 
into slow modes, |q| < go ^ Pf- (a), (b), and (c) are the 
Hartree, tlie Fock, and the Cooper vertices, respectively. In 
our model, the Cooper vertex should be omitted to avoid 
double-counting. 



modes of the microscopic model, Eqs. (|2.4p - (|2.6|) . Effec- 
tively, only those fermions contribute that are energeti- 
cally close to the Fermi energy ep- From this constraint, 
wc obtain the three relevant vertices shown in Fig. [1] de- 
scribing scattering processes with momenta p, p' located 
near the Fermi surface. We single out these vertices as- 
suming that the momenta q are small, |q| < goi where go 
is a phenomenological momentum cutoff which is much 
smaller than the Fermi momentum pp . Vertices (a) and 
(b) describe soft interactions in the particle-hole chan- 
nel, whereas (c) is the Cooper vertex. The well-known 
Hartrce-Fock approximation is obtained using the ver- 
tices (a) and (b) and, thus, they can be referred to as 
the Hartree and Fock vertex, respectively. The Cooper 
vertex enters the ladder diagrams leading to the BCS 
superconducting instability in case of attractive interac- 
tion. 

In principle, all three vertices are important when 
calculating physical quantities. However, we are here 
interested in the anomalous contributions to the ther- 
modynamics, which originatei^iHii^ from (quasi-)one- 
dimensional processes. As we will see later, the main 
contribution to the anomalous terms comes from small 
q, which is also seen from the conventional perturba- 
tion theoryfiS Of course, assuming that all the vertices 
in Fig. [T] arc different forbids the regions of the essential 
momenta attributed to them to overlap. 

A quick glance, however, reveals that this is not the 
case. For example, there is a clear overlap between the 
regions of the momenta in the vertices (a) and (b) at 
small |p — p'l ^ |q|. In order to avoid double-counting, 
one would have to consider only one of these vertices in 
this region. In Ref. [l^, e.g., it was chosen to remove the 
region |p — p'| ^ |q| from the Fock vertex. Fortunately, 
the contribution coming from this region of the momenta 
can be neglected at the low temperatures considered here. 
As a result, we can consider the vertices Fig.ma) and (b) 
as practically different vertices for |q| < go ^ Pf- 

At the same time, the Hartree-Fock (a,b) and the 
Cooper (c) vertices do overlap in important momentum 
regions. As to the role of the Cooper vertex, building 
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ladders out of it, logarithmic divergencies appear for ar- 
bitrary scattering angles pp'. The contributions coming 
from large angles are not reproduced by using the Fock 
vertices (b) instead. However, scattering angles essen- 
tially different from or tt are less important for the 
anomalous contributions we wish to calculate. Consid- 
ering only angles pp' close to or tt means focusing on 
almost one-dimensional scattering and this is where the 
anomalous contributions emerge. 

This region of the momenta attributed to the Cooper 
vertex Fig. [IJc), however, fully overlaps with that for 
the Hartree-Fock vertices (a,b) . Any diagram containing 
Cooper loops with small angles pp' can be represented 
in an equivalent way using particle-hole loops built from 
vertices (a) and (b) . Several examples of this equivalence 
can be found in Ref. [isl Therefore, taking into account 
all the vertices Fig. [Ha,b) and (c) would imply double- 
counting. In order to avoid it, one should choose between 
either the Hartree-Fock vertices or the Cooper ones, but 
not take into account all of them. 

In the present study, we choose as in Ref. [l^ the 
Hartree-Fock route. This is in contrast to the approach 
of Ref. [l3 where the Cooper channel representation was 
used. Our choice will turn out more convenient for sin- 
gling out the anomalous contributions. 

As a result, we write the effective interaction describing 
the low energy physics as 

'5int = ^ J2 {x*AP + Q)xAP)V^X*a'{P')Xa'iP' + Q) 
PP'Q,aa' 

- xUP + Q)Xa-{P)V^-p,x'AP')Xa{P' + Q)] 

(2.7) 

instead of Eq. (|2.6p . In Eq. (|2.7p . the fermionic fields are 
represented in Fourier space and four-momentum nota- 
tions 

P=(e,p), Q = (w,q) 

are used. Herein, e and w are fermionic and bosonic 
Matsubara frequencies, respectively. In the "fermionic" 
summations Ep(')(- • ■) = T^EeO /(■ • •)[d''P^'V(27r)'^], it 
is understood that p, p' arc of order pp- On the con- 
trary in the "bosonic" summation X)q' introduce a 
function /(q) which cuts off the momentum q beyond 
go < Pf, 

E(-)-J^E/(-)/(q)(g,- (2.8) 

The function /(q) can for instance be modeled as /(q) = 
0(90 — |q|) with 6 denoting the Heaviside function. In 
final formulas, we may choose a different form which al- 
lows to conveniently perform the terminal integrations. 
For the moment, we do not specify it more than that it is 
assumed to fulfill /(q = 0) = 1 and decay fast beyond ^q. 



C. Bosonization 

Our choice of the effective interaction, Eq. (|2.7p . leads 
after bosonization to a field theory representing particle- 
hole-type bosonic excitations of the Fermi gas. The route 
to obtain the effective low energy field theory for these 
excitations follows the higher-dimensional bosonization 
scheme introduced in Refs. and (2^ It is done in five 
steps: 

1. Decouple the effective interaction employing a 
Hubbard-Stratonovich transformation. 

2. Integrate out the fermionic degrees of freedom. 

3. Derive the effective equation of motion for the 
bosonic field subject to the random Hubbard- 
Stratonovich auxiliary field. 

4. Write the solution of this equation in form of a 
functional integral over superfields thus obtaining 
a closed supersymmctric field theoryi^"— 

5. Average over the auxiliary field. 

In what follows, each of the steps is presented in a sepa- 
rate section. 



1. Hubbard-Stratonovich transformation 

In order to decouple the effective interaction iSint by 
an integration over the auxiliary Hubbard-Stratonovich 
field, wc recast the spin structure in the Fock channel by 
means of the relation 

3 

where (ct^) = (o-^,"") with = 1 and cr = (o-\(T^,cr^) 
is the vector of Pauli matrices. In what follows, we use 
the convention that Greek upper indices appearing twice 
imply the summation from to 3. 
Wc thus obtain 

Q 

-^Y.^P-P'S^iQK'{-Q)} (2.10) 

pp' 

where 5*^ are components of the four-component vector 
{Sp, Sp, Sp, Sp) and 

n(Q) = E^-(^'P + f)^-(=-^'P-f)' (2-11) 

P,(T 

S^iQ) = T ^ x:(£, P + ^K.'X.' (£ - ^, P - |). 
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The zero component is related to the particle density 
n{Q) as 

Y,S°{Q) = n{Q) (2.12) 
p 

whereas the other components form the spin density vec- 
tor S(Q), 

^Sp(Q) = S(Q) (2.13) 
p 

Assuming that the interaction V (r) decays sufficiently 
fast, the interaction amplitudes V entering Eq. (|2.10p can 
be written as 

V^^Vo, Vp_p. ^y(2pf^sin(^)) (2.14) 

with pp' denoting the angle between p and p'. In 
Eq. (|2.14p . Vq is a positive constant and F is a func- 
tion of the angle between the vectors p and p' located at 
the Fermi surface. As the main contribution will come 



from momenta p near the Fermi surface and small q, 
Eq. ()2.14|) is a good approximation. For contact interac- 
tion, V = Vo. 

In principle, we are now ready to decouple the quar- 
tic interaction, Eq. (|2.10p . by means of a Hubbard- 
Stratonovich (HS) transformation. Considering the low 
energy theory, we have not said anything about high en- 
ergies so far. Actually, one can first integrate out the 
high energy degrees of freedom. Following the philoso- 
phy of the Landau Fermi liquid theory, integrating out 
the high energies results in a renormalization of the origi- 
nal spectrum, Eq. (|2.2p . of the interaction, Eq. (|2.3p . and 
of the ground state energy. Neglecting fluctuations, one 
would obtain a partition function containing the renor- 
malized constants to be considered as phenomcnological 
parameters. 

As a result of the integration over the high energies, 
one obtains instead of the quadratic forms nn and S^S^^ 
written in Eq. (|2.10p expressions of the type (n— (n))(n — 
(n)) and (5'^ - {S^'))iS^' - {S^)), where (...) denotes 
the averages with respect to the high energy part of the 
Hamiltonian. Applying the HS transformation to this 
effective interaction yields 



exp{-5i„t} = / exp{i^^(-Q)[n(Q)-(n(Q))] +^^/i^;(-g)[5^(Q)-(5(i(Q))]} W^[^]W/[/inP^P/i'' . 

Q Q P 

(2.15) 

I 



The Gaussian weights entering Eq. ()2.15p are given by 

=exp{-^^V'(Q)^(-Q)}, (2.16) 
Q 

Q pp' 

and VipDh^^ is the measure normalized such that 
/ >VftM>V/[/i'']X>(pm^ ^ 1. The HS fields ip{Q) 
and h'^iQ) depend on bosonic Matsubara frequencies w, 
corresponding in imaginary time representation to peri- 
odicity, (p{T+l3; q) = ip{T; q) and /i^(t-1-/3; q) = hi^j r; q). 
The constant Vq has been introduced in Eq. (|2.14p and 
is the inverse of the operator Vp which acts on a 
function g{p) as 

[K9](P)-EW.9(P')- 
p' 

Actually, V^-p' is essentially the function V from 
Eq. (|2.14p for momenta p close to the Fermi surface. By 
construction, the momentum q in the HS fields (p{Q) 
and h!^{Q) is restricted by the cutoff qo Pf and the 
coupling of the HS fields to the fermions displaces the 
fermion momentum p of order pp only locally. 



2. Integration over the fermionic fields 

After the HS transformation, Eq. (|2.15p . the partition 
function is given by 

Z = Zq j W[$]2?<I> (2.17) 

where $ is the 2 x 2-matrix 

*p(o) = i^(g) +/l^^(Q)'T^ (2.18) 

while W[$] WhMWf[h^'] and I>$ = VipVht'. The 
functional Z[^] is essentially the formal partition func- 
tion of non-interacting fermions in the field 3>, 

Z{^] = Z^'exp{-T[<^>]} J exp{ -cSoi*]} V{x*,x) 

(2.19) 

with the action <So [$] given by 

^50 m = 5o - r 5] xU^. P - f ) P q)x (r, P + f ) dr. 

pq 

(2.20) 
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Herein, the bare action So has been defined in Eq. (|2.5p , 
the field $ is written in the imaginary time r representa- 
tion, and the spinor notation x ~ (xt'Xj.) is used. The 
functional J-"[$] in Eq. (|2.19p is given by 

m= rE(x^(^'P-f)*p(^'q)x(^.p+f))dr 

" pq 

(2.21) 



Herein, is the Green's function of the ideal Fermi gas 
described by the action Sq , Eq. (12.51) , and tr denotes the 
trace over spins. 

In order to integrate out the fermion fields in 
Eq. (j2.19p . we follow the route suggested in Refs. Hm 
and recast the fermion determinant as det = expTrln. 
The logarithm is then replaced by an inverse using an 
additional integration over a variable u of the form 



ln(a; + (/))= In X + / (pix + uq 
'o 



This yields 



(2.22) 



=exp{-J'[$]} 



where is the 2 x 2-matrix Green's function of non- 

interacting fermions in the external field u$, i.e. the 
propagator of the action iSo[u$], Eq. ()2.20p . It satisfies 
the equation 



= S{t - r')5q,o . 



2$] 



(2.23) 



Thus, in order to study the thermodynamics of the low- 
lying excitations, we need to determine the Green's func- 
tion of non-interacting fermions in a random external HS 
field at equal times. This quantity (r, t+O) is a com- 
plex matrix function being periodic in imaginary time r. 



a["*l(T-|-/3,r-|-/3 + 0) = a["*l(T,-r-hO). This observation 
motivates the introduction of the complex 2 x 2-matrix 

field ^pp- = A^;^} as 



-4pp'(r) 



g^l{r,r + 0)-g^:':Hr,r + 0). (2.24) 



The field ^pp'(t) obeys bosonic periodicity in time r, 
-4pp' {t + l3) = -4pp' (r) , and captures the entire physics 
of the low-lying excitations. 

By this definition, the functional Eq. (P?^ . 

takes the form of a functional of the boson field .4pp'(r), 
Eq. (EUl), 

Z['f]=cxp{- [ r^tr[<i>p(T,q)^j^.,j^. (r)] drdu} . 

^ -^OJO pq J 

(2.25) 

Using Eqs. (|2.17p and p.25p one can compute the parti- 
tion function Z provided the function ^pp' (r) is known 
for any configuration of $. Of course, one could solve 
Eq. (|2.23p for the Green's function and find -4pp' (r) 
from Eq. p.24p but this would lead to the conventional 
diagrammatic expansion for the fermions. Here, we fol- 
low a different route deriving a closed equation for the 
field ^pp- (r). 



3. Dynamics of excitations 



The derivation of the equation of motion for the boson 
field Apr)' follows the route presented in detail in Refs.[23l 
and[2J: We start with the differential equations for the 
Green's functions g'^^^T^r') and g^'"^\T,T'), Eq. (I^?^ . 
Subtracting them according to the definition of ^pp', 
Eq. (|2.24p . yields a differential equation as an interme- 
diate result. Next, we repeat the same steps with the 
equations for ^/["'(t, r') and ^^"*i(r, r') while this time 
the operator dr + Ho — acts from the right. The 
equation obtained in this way is then subtracted from 
the intermediate result. Eventually putting t' = t + 0, 
we find that the dynamics of .4p^k p_k is governed by the 
equation 



dr 



(p.k)- 



A 



pft,P-| 



I>ff4<J,p-f 



-4pf|.p-§-<j*I^Ma(q) 



-k)[r 



(2.26) 



Herein, Up ~ [exp(^p/r) + 1] ^ with ^p ~ e{p) — is 
the Fermi distribution function of the ideal Fermi gas. 

In Eq. p.26p . the relevant momenta k in Ap^k p_k 
arc small and do not exceed the cutoff momentum qo <C 



Pf- This can easily be seen from Eq. (|2.25p because the 
fields $p (t, k) arc by construction non-zero only if the 
momentum k is small. The fact that the two momenta 
p + k/2 and p — k/2 entering -4p^_k p_k are essentially 
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close to each other significantly simplifies the analytical 
study as compared to the general formulation studied 
perturbatively in Ref. [13 as a check of the bosonization 
procedure. 

In order to obtain a feasible low energy theory, the 
still present degrees of freedom on the scale pp should 
be integrated out from Eqs. (|2.25|) and (j2.26p . Since the 
right-hand side of Eq. (|2.26p contains the combination 
n„ k — n„ I k . the dependence of the function ^^k ^k 
on IpI near pf is for small k sharper than the depen- 
dence of all other functions entering the left-hand side 
of Eq. (|2.26p or the exponent in Eq. (|2.25p . As a result, 
when integrating both sides of Eq. (|2.26|) or the integrand 
in the exponent of Eq. (|2.25p over the variable ^p, it is 
justified to approximate p ~ pfh with = 1 in all other 
functions depending smoothly on |p|. 

Proceeding in this way, we rewrite the exponent in the 
expression for Eq. (j2.25p . as 

5]$p(q)^j^a,p_s (2.27) 
pq 

thus separating the radial and angular integration of the 
momentum vector p = [2m(ei? -I- ^p)]^/'^n. In Eq. (|2.27p . 
V is the density of states at the Fermi surface and the 
integration /dn is done over the (d — l)-dimensional 
sphere S''"^. We normalize the integration over n by 
the convention 



dn — 1, 



(2.28) 



and write the function $p (q) as $„ (q) since p pr'n.. 

Equation (|2.27p shows us that we need the integral 
over fo of ^^k ^k rather than this function itself. This 

^ Pi 2 'P^2 

observation motivates us to introduce the quasiclassical 
field 



an(k) 



A 



pff.p-4 



d^p 



and to construct the low energy theory for this field. 

For this purpose, we integrate the equation of motion 
for ^p^k p_k, Eq. (|2.26p . in a similar manner term by 

term. For the first term of Eq. (|2.26p we find 



dr 



(p.k) 



A. 



pff .p-^ 



;d^p ~ 9r + VF{n-'k) an(k), 
(2.30) 

where is the Fermi velocity. 

For the first expression in the second term, we obtain 



pfii±a+a,p_is+a+^d^P 



(2.31) 



^^iLLiSi {(l)a-n+^ (k + q) 



2pp " ' 2pp 

where k^ is the component of the vector k perpendicular 
to the vector n. 



k (n • k)n 



(2.32) 



In Eq. (|2.3ip . we made the same approximation as in 
Eq. (j2.27p using the fact that the field $ is a slow function 
of the variable |p|. The radial integration in Eq. (|2.3ip 
absorbs the normal component of q and that is why only 
the tangent component q_L enters the result of the inte- 
gration. The combination n + q±/{2pp) in Eq. (|2.3ip 
for small vectors q^ corresponds to a rotation of n. Ac- 
counting for these rotations is very important because 
they capture the essential effects of the curvature of the 
Fermi surface arising in d > 1. 

The remaining terms of Eq. (|2.26p are treated simi- 
larly. Integrating the right-hand side is especially simple 
because ?T.p_k — ^p^k ~ m~^(p-k)(5(^p). Finally, we re- 
duce the functional Eq. (|2.17p . to the form 

= exp{ - y" ^tr[$„(q)a„(q)] dndrdu} , 

(2.33) 



(2.29) where an(k) is the solution of the equation 



dr+VFin-k) a„(k) - tty^|c& k^-w^ (q)Q^+q^ (k-Hq) -a„_ (k+q)$ k^^^ (q) = -ui/t;F(n-k)$„(-k) 



(2.34) 



Let us discuss the analytical properties of the 
field an(k) = an(M,T, k) and its equation of motion 
Eq. (|2.34p . By construction, the field an(T, k) describ- 
ing the low energy excitations is bosonic: It is a complex 
2 X 2-matrix field periodic in time, an(T-|-/3,k) = an(T,k). 
The momentum k entering the field a, |k| ^ 90 ^ Pf, 



determines the scale of the bosonic excitations while the 
argument n of the field a„(k) determines the position on 
the Fermi surface. The dependence of an(k) on n and k 
contains the full information needed to describe the spin 
and charge excitations in a higher-dimensional Fermi gas. 
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In dimension d = 1, transverse momenta q_L, do 
not exist. If we neglected these momenta in higher di- 
mensions, we would come to the low energy model of 
Ref. [l^. In this approximation, Eq. (|2.34p describes one- 
dimensional processes and all effects of the Fermi surface 
curvature in d > 1 are lost. Writing the field an (k) in 
the form an(k) — £in(k)cr'' + Sn(k) • cr, one can sec that, 
as a result of this approximation, the charge Qrx (k) and 
spin Sn (k) parts of the field decouple from each other. 
Eventually, one obtains a model of non-interacting charge 
excitations gk(n) and a non-trivial field theory for the in- 
teracting spin excitations Sn(k). All the results of Ref. [isl 
have been obtained in this way, implying their validity in 
one dimension but requesting reconsideration for d > 1. 

Taking into account the curvature effects in dimensions 
d > 1, Eq. (|2.34p shows that the spin and charge excita- 
tions cannot be treated separately but interact with each 
other. Moreover, the charge modes also interact them- 
selves similarly to the spin modes. The interaction of 
spin and charge modes is what constitutes the major dif- 
ference between the physics of the fcrmion gases in d = 1 
and d > 1. Furthermore, we will see that there are classes 
of diagrams in a perturbative analysis of the final boson 
theory that are logarithmic in d = 1, but due to the 
presence of the q_L-terms become regular in dimensions 
d > 1. 

In principle, one can solve Eq. ()2.34|) employing a per- 
turbation theory in the HS field $. In the zero-order 
approximation, one should neglect the terms containing 
the HS field $ in the left-hand side of Eq. (HSJ), yielding 



a„(w,k) 



,(0) 



(w, k) with 



iw — Wi?(n-kj 



(2.35) 



Inserting this zero order approximation an''('^,k) into 
the functional Eq. (|2.33p . we reduce the partition 

function Z, Eq. (|2.17p . to a Gaussian integral over the 
field $ = \ip+h^a^, Eq. (EHH]). It is not difficult to under- 
stand that this limit yields the contributions obtained by 
summing certain ladder series in the conventional fermion 
diagrammaticsj^ Considering only the HS field while 
neglecting h}\ we obtain the contribution of the rings 
built from polarization bubbles, i.e. reproduce the ran- 
dom phase approximation (RPA). Alternatively, keeping 
only W the contribution of the particle-hole ladder ring 
is reproduced. Keeping both and /i^, we obtain the 
contribution of all particle-hole ring diagrams. 

However, the interesting logarithmic contributions to 
the non-analytic terms arise from the fluctuations on top 
of these ladders. Considering the ladders as propagators 
of elementary bosonic excitations, one can say that the 
logarithmic contributions arise as a result of interaction 
between these excitations. In order to treat the fluctua- 
tions properly, we need a more efficient route of solving 
the equation of motion for an(t^, k), Eq. (|2.34l) . including 
the $-term in the left-hand side. 



4. Superfield representation 

In this section, we represent the solution of the equa- 
tion of motion for the bosonic field a„(k), Eq. (j2.34p . for 
an arbitrary $ in the form of a functional integral over 
superfields. We begin by noting a remarkable symmetry 
in the left-hand side of Eq. 

Being linear in an(k), the left-hand side can be for- 
mally represented as [£a]n(k). We observe that the oper- 
ator C is antisymmetric with respect to the inner product 
given by 



(/^' 



= /WEt^[/n(-k)an(k)] dndrdii. (2.36) 
Jo Jo J 



Herein, is a field having the same structure like a. The 
antisymmetry condition {f \ Co) — —(a, is straight- 
forwardly checked using the definitions of the operator C 
and the inner product Eq. (|2.36p . We will see shortly 
that this antisymmetry of C leads to an important sim- 
plification of the theory. 

Since the remaining of the derivation of the low en- 
ergy field theory is purely formal, we will use short-hand 
notations in this section. We define 



n = n^iu, T, k) = -uzyi>F(n-k)$„(-k) 



(2.37) 



as short-hand notation for the right-hand side of 
Eq. (|2.34p and use the notation C introduced above for 
the antisymmetric operator in the left-hand side. As a 
result, Eq. (|2.34p can be written in the compact form 



(2.38) 



Both the operator C and the inhomogeneity term TZ de- 
pend (linearly) on the HS field Eq. dUg) . W ith 
$ being a Gauss-distributed random field, Eq. (|2.38p is 
technically a stochastic differential equation for the bo- 
son field an(w, r, k). In the context of stochastic field 
equations, a well-known method of analysis is the Bccchi- 
Rouct-Stora-Tyutin (BRST) transformation^^"— which 
brings the problem of solving the stochastic equation into 
the form of a supersymmetric field theory. The latter for- 
mulation allows for a study by means of standard field 
theory techniques. 

We now apply the BRST map on our problem. First, 
we rewrite Eq. (|2.33p with a satisfying Eq. (|2.34p in a 
form of a functional integral over fields a 

5\Ca - n] dct— Z\-^-a]Va^Va. (2.39) 
da 

Herein, the integration with respect to the measure 
Da'^'Da is performed over all complex fields a which 
do not necessarily satisfy Eq. (|2.34p . The functional 
Z[<I>; a] in Eq. (|2.39p denotes formally the functional from 
Eq. ()2.33p . yet the field a is here an unspecified com- 
plex field and not the solution to the constraint equa- 
tion Eq. (j2.34p . The equivalence of the representation by 
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Eq. ((OQ)) and the original one, Eqs. ((2:33)) and (|234)) . for 
Z [$] is easily seen as the constraint (|2.34p is enforced by 
the integration over the functional 5-function in the inte- 
grand of Eq. (|2.39|) . The determinant in Eq. (|2.39p arises 
as a consequence of changing variables from a to Ha. 

Our goal is to integrate the functional Z [$] , Eq. p.39p , 
over the fields $ and obtain a field theory for the inter- 
acting bosonic excitations. This can be achieved repre- 
senting the (5-function as a Fourier integral, 

8\La - n] (2.40) 

exp {'-{f\Ca~n)+^{[Ca-n]\.f)] VfVf , 

and the determinant as an integral over Grassmann vari- 
ables. 



det 



SC 
6a 



(2.41) 



oc 



exp{^ {p\Ccj) + ^ {[Ccj]\p) } Va^VaVp^Vp 



In Eqs. (|2.40p and (|2.4ip . / is a complex field, while a 
and p are Grassmann fields of the same structure as a. 
The brackets (...,...) denote the inner product defined 
in Eq. ^M^. All the fields in Eqs. and (f2:4T|) 

are assumed to be periodic in imaginary time r in or- 
der to reproduce the bosonic boundary condition of the 
solution a of Eq. (|2.34p . 

Substituting Eqs. (PTiO)) and (piiT)) into Eq. 
we come to the representation of the functional 
Z [$] in terms of a Gaussian integral over the fields 
a,a\f ,f\a,(j^ ,p,p^ . However, as the functional Z[$;a] 
contains only the field a and not a^ , the integral over 
the fields a\ f,a\ p can immediately be calculated giv- 
ing unity. Then, we are left with an integral only over 
the fields a, , cr, . 

Instead of writing all these fields separately we unify 
them into one superfield 5* which we define as 



*( 



(2.42) 



9 and 6* are additional Grassmann anticommuting vari- 
ables. By construction, ^' is an anticommuting field. 
This, however, does not mean that it describes fermions 
as the periodicity in imaginary time, 



4- (r) * (r -f /?) 



(2.43) 



guarantees the boson statistics. 

The antisymmetry of the operator £, Eq. (|2.38p . im- 
plies the remarkable and important relation 



j («-, £*) Aede* = 2 {f\Ca) + 2 {p\Ca) 



Using additionally the relations 



(2.44) 



(2.45) 



we can express the entire field theory solely in terms of 
the superfield 'J. 

As a result, we write the partition function Z of the 
low energy field theory for the excitation modes of the 
interacting Fermi gas, Eq. (j2.17p . in the form of a func- 
tional integral over the superfield ^' and the auxiliary 
field 

Z = Zo J exp { - 5s - 5b } W[^]V<^>V'if . (2.46) 



The action Ss with 



Ss 



{'^,C^) + i'^,'R)0 dOdd* (2.47) 



is invariant under the BRST symmetry transformation 
^' M> ^' + (5\1/ with the variation given by 



9^ 



{2.4 



Herein, the transformation parameter 77 is a Grassmann 
variable. The action Sb, which derives from the func- 
tional a], takes the form 



1 rP 



Jo 



^tr[$n(q)*n(q)]6'* dSdrdndrdM 



(2.49) 

In contrast to Ss, Eq. (j2.47p . the action Sb is not invari- 
ant when varying the superfield ^' according to Eq. ()2.48p 
and, thus, breaks the BRST symmetry!^ 



5. Final form of the low energy model 

Since the action Ss + Sb is linear in the auxiliary 
field <&, the integral over $ in Eq. (|2.46p is Gaussian and 
can easily be performed. This yields the final form of 
the low energy field theory for the bosonic excitations of 
the interacting Fermi gas. All the interesting physics is 
described by the 2 x 2-matrix superfield \1/ only. The par- 
tition function can be written after the integration over 
the HS field $ in the form 

Z = Zo y exp { - 5barc - 52 - 5,3 - 54 } . (2.50) 

Herein, the bare action 5barc is in Fourier representation 
given by 



5b 



i /'^tr[*„(-/^,«;) (2.51) 

K 

{-ie-|-i;F(n-k)}^'„(A', k)] dndK . 
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For the sake of compact notations, we use here and in 
the following the four-momentum notations 

A' = (s,k),^(...) = r^/(...)|^, 



Sk.-k' 



T 



k,-k' 



Q=.(.,q),B-)-^E/(-)7^, 

where e and w denote bosonic Matsubara frequencies. 
Also, we use the following short-hand notations in the 
remaining of our analysis: 



K = iu,e,e*) , 

dn = d9de*du , 

s{k - k') = 6{u - u'){e* - e'*){d - 9') 



(2.52) 



Whenever we integrate over u, integration over the inter- 
val (0, 1) is implied. 

Averaging quadratic forms with respect to the bare 
action 5baro7 Eq. (|2.50p . is done as 



(2.53) 



x5„a'Saa'S{K. - K')6{n-n')5K-K' ■ 
In Eq. (1^1) . 

9n{K) ^ 



ie — vp{n-'k) 



(2.54) 



is the bare Green's function for the bosonic modes. 
Higher moments of the field 4* are reduced to second 
moments, Eq. (|2.53p . using Wick's theorem. 

Let us now have a look at the interaction vertices in 
Eq. p.50p . The quartic interaction term 54 is given by 

^4 = J dndndftdK ^ uu 7j^/(q) (2.55) 



KKQ 



X tr 
X tr 



*fi(-/?,K)vI/^_^(/?-Q, 



^ j dndndKdK E uu 7|^/(q) 



KKQ 



X tr 
X tr 



*„(-i^,«:)a^-*„+ji^(A' + Q,K) 



*ii(-X,K)a**._j^(A'- 

2PF 



where the vector is the projection of q onto the plane 
perpendicular to n, i.e. q^ = q— n(n-q). The amplitudes 
for the spin and charge channel are expressed at 

nh nh 

weak interaction in terms of the interaction potential V ^ 
Eq. (EUl), as 

7:-. = -in2pi^-n^)' 7:~, = ^^o+7„^,, (2.56) 



respectively. If the interaction is not weak these ampli- 
tudes can be considered as effective coupling constants 
of the Fermi liquid. The cutoff function /(q) introduced 
for the soft modes in Eq. (|2.8p is from now on written ex- 
plicitly in the formulas. The cubic interaction reads 

S3 = j^f dndndAtdS^ 7. 7^/(q) (2.57) 

KQ 

X tr *„(-A^k)^'„+jll.(A' + Q,k) 

X ^^~uvvF{n■q)e + 210* | tr [*fi(-g, k)] 

+ dndndKdK^u 7^/(q) 



KQ 



X tr 



X i^-uiyvF{n-q)e + 2i0*^tT [cr''*fi(-Q, «)] 

and, finally, the quadratic action 52 takes the form 

" /dndfidAcd^^ 7£./(q) 

{uiyvF{n-q)9 + 2ir } tr [*„(Q, k)] 
|-wwj.(n-q)6i + 2ir } tr [*s(-(9, k)] 

i- / dndridAcd;^^ 7^/(q) 

{uiyvF{n-q)9 + 2ir } tr [a'^^'niQ, n)] 
\^-uiyvF{n-q)0 + 2ir | tr [cr''*fi(-Q, k)] 



(2.58) 



X 

X <; - 



X 

X <; - 



Diagrammatically, pcrturbative calculations within the 
low energy boson model can be conveniently represented 
using the building blocks shown in Fig. [2l In the di- 
agrammatic representation and explicitly in Eqs. ()2.55p 
and (|2.58p . it is evident that the building blocks consti- 
tuted by 1S4 and ^2 are invariant under vertical reflection. 
On the other hand, it is important to note that for d > 1, 
neither ^4 nor S3 arc symmetric with respect to horizon- 
tal hipping. 

The quartic interaction 1S4, Eq. (|2.55p . is fully ob- 
tained from the ^-average of the BRST-symmetric ac- 
tion Ss, Eq. (PTf)) . Since ^4 inherits this BRST sym- 
metry, there is no (perturbative) contribution to thermo- 
dynamics originating purely from ^4. In contrast, the 
cubic and quadratic interactions 1S3 and ^2, Eqs. (j2.57p 
and (|2.58p . are formed using also the symmetry breaking 
action Sb, Eq. (|2.49p . when averaging over the auxiliary 
field $. That is why diagrams contributing to a phys- 
ical thermodynamic quantity necessarily contain ^2 or 
1S3 among their building blocks. For example, consid- 
ering only the terms 5barc + ^2 we reproduce the RPA 
and particle-hole ladder rings because this is equivalent 
to neglecting the auxiliary field $ in the left-hand side 
of Eq. (|2.34p . cf. the discussion at the end of Sec. Ill C 31 
Blocks of 54 may additionally decorate diagrams built 
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(a) gn{K) 



-if,n 



(a) 



(b) 



54 



K-Q,h- 



-K,n 



2p,.- 



qi 



2pp 



FIG. 2: Diagrammatic building blocks of our low energy field 
theory: (a) the propagator (?n(A"), Eq. (|2.54[l . and (b) the 
interaction vertices <S4, Ss,, and 52 from Eqs. (|2.55|l . (|2.57fl . 
and (f238|l . 



from 52 or and the contribution may consequently 
acquire logarithmic renormalizations. 

In summary, Eqs. (|2.50p - (|2.58p specify our effective 
low energy model for the interacting Fermi gas in d > 1 
dimensions. It is a field theory for the anticommuting su- 
perfield 4* which describes the bosonic excitations. The 
interaction between these excitations appears as sum of 
the quadratic term 52, the cubic term ^3 and the quar- 
tic term ^4. The bare coupling constants are written 
in Eq. (|2.56|) . In principle, one can immediately start 
perturbative studies of the model using the contraction 
rule, Eq. (|2.53p . and Wick's theorem. A possible dia- 
grammatic representation is shown in Fig. [2] Although 
the effective field theory may look somewhat complex, it 
allows to conveniently treat the low energy limit, identi- 
fying the interesting logarithms and summing them. This 
is what the next sections are devoted to. 



Q,n 



(b) 



54 



(c) 



ni 



"2 



FIG. 3: Backscattering diagrams for the thermodynamic po- 
tential Q.: Diagram (a) is the second order diagram containing 
the leading backscattering contribution for n ~ —ii while di- 
agram (b) represents an exemplary logarithmic renormaliza- 
tion to diagram (a), cf. Sec. IIII Cl Finally, diagram (c) gives 
for ni ~ —j\2 ~ — ng ~ n4 a backscattering contribution 
of higher order in the interaction and also includes a renor- 
malizing building block <S4. For weakly interacting fermions, 
diagram (c) can be neglected. 



in d = 1 dimension are suppressed in dimensions d > 1 
because of transverse fluctuations c\± along the Fermi 
surface. Eventually, the effects of the finite Fermi sur- 
face curvature lead to renormalization group equations 
different from the ones obtainedi^ in one dimension. 

Before studying the one-loop vertex corrections, we be- 
gin the perturbative analysis of this section as we dis- 
cuss the relevant diagrams for the thermodynamic po- 
tential. These diagrams describe physical backscattering 
•processes. 

While the boson model, Eqs. ([^30)) - ([^35)1 . has been 
derived for an arbitrary dimension d, we consider from 
now on the most interesting case of a two-dimensional 
Fermi liquid, d = 2. 



III. PERTURBATION THEORY 

The bosonized model, Eqs. (|2.50p - (|2.58p . is not triv- 
ial and the perturbation theory in the coupling constants 
75-:, 7^:;, Eq. (j2.56|) . yields logarithmic contributions di- 
verging in the limit T — > 0. In this section, we iden- 
tify the relevant classes of logarithmic one-loop diagrams. 
Later in Sec. IIVI these logarithmic contributions will be 
summed up to infinite order by means of a one-loop renor- 
malization group scheme. 

In one dimension, such a procedure would essentially 
repeat the steps from Ref. [Isl. The peculiarity of higher 
dimensions, d > 1, appears in form of the "rotations" 
n (\ji_/2pp of the angular arguments in the interact- 
ing superfields, cf. Eqs. (|2.55p - (|2.58p . Consequently, the 
running momentum Q in a one-loop diagram affects at 
the same time the (actual) momentum K and the di- 
rection n of the propagators. As a result, we will find 
that logarithms which certain classes of diagrams feature 



A. Backscattering diagrams 

In the second order in the interaction, only diagram 
Fig. [3l^a) describes a contribution to the thermodynamic 
potential relevant for studying the backscattering ef- 
fects. All other second order diagrams cannot contain 
two boson propagators gii{K) and gn{K') with n ^ — n. 
Figure [S^b) shows an exemplary diagram that renormal- 
izes the bare diagram Fig. [3ja) while Fig.[3][c) represents 
a backscattering contribution of higher order in the in- 
teraction. Considering the limit of weak interaction, we 
are safe to neglect such higher order diagrams because 
they do only describe high energy renormalizations of 
the coupling constants. 

Working with the effective low energy theory, wc have 
to be sure that the main contribution to the physical 
quantities of interest indeed comes from the low energies 
not exceeding T. Whether this is the case or not, it 
should be checked for each quantity under investigation. 
In fact, the low energy contributions arc not the most 



12 



important for a perturbative correction Afi (T) to the 
thermodynamic potential and, thus, we cannot compute 
Aft (T) using the low energy limit only. However, the 
main contribution to the difference 



sn{T) = An{T) - An{T = o) 



(3.1) 



does come from the low energies. In order to determine 
such quantities as the specific heat, the quantity Sfl (T) 
contains all the necessary information and the low energy 
bosonized model, Eqs. p.50|) - (|2.58p . becomes useful. 

Formally, the quantity AD, (T) will be represented 
in terms of sums over Matsubara frequencies such as 
TJ2uj '0('^n)- The corresponding expression for Sn{T) 
consequently takes the form 



2^ 



E 



-iujl/T 



dw 
2^ 



(3.2) 

Equation (|3.2p follows from the Poisson summation for- 
mula. It shows that, when calculating dVl (T), essential w 
in the function ip (w) are of order T provided the function 
tp (w) decays sufficiently at \uj\ — > c». 

Using the developed formalism, we can start the cal- 
culation of thermodynamic quantities. As a first exam- 
ple, we are going to compute the correction JSl'^' (T) in 
the second order in the coupling constants 7?- and 7S- . 

nn nn 

Logarithmic contributions arc taken into account later by 
replacing the bare coupling constants 7^ and 7^ with 

nn nn 

effective amplitudes obtained from the summation of log- 
arithmic contributions. This computational procedure is 
justified by the fact that logarithmic contributions come 
from energies exceeding T with the logarithms cut from 
below by max{2TrT,vpq()\n + n|). Therefore with loga- 
rithmic accuracy, one may replace the energies in the ef- 
fective amplitudes by the temperature T and treat them 
as constants when calculating Matsubara sums. 

The second order contribution is given by the dia- 
gram Fig. [2Ia) . Analytically, we obtain for An^^^{T) = 
— (T/2) the expression 



correction (50^^^ (T) as 




(27r)2 
VF{n-q) vp{h-q) 



dndn /2(q) 



X7^-; 

[io; — vp{n-q)][—iuj + vp{h-q)] 



(3.4) 



The scattering amplitudes for the charge and spin chan- 
nels, j'^ and ji-: , enter Eq. p.4p as 



' nn^ 



nfi^ 



(3.5) 



The evaluation of the integral in Eq. p.4p follows with 
some minor deviations the steps of a similar calculation 
in Ref. [isl . This calculation is not trivial and we present 
a possible route of how to carry it out in Appendix [X] 
As a result, we find 



TTVr 



(3.6) 



By virtue of the thermodynamic relation Sc = 
—Td^5il/dT^, the anomalous correction to the specific 
heat 6c is in the second order in the interaction obtained 



as 



(3.7) 



We see that only the backscattering amplitude [nn = tt] 
enters Sc'^'^^T), Eq. (ir7| . 

Equation p.7p gives the well-known anomalous lowest 
order specific heat contribution.— It is quadratic in T, 
which contrasts what one would expect from the Som- 
merfeld expansion for the Fermi gas of weakly-interacting 
quasiparticles. Thus, Eq. p.7p confirms the equivalence 
of the perturbative calculations in both the conventional 
approach and the bosonization one which we are study- 
ing here. In the remaining of this paper, we investigate 
the logarithmic rcnormalizations to the backscattering 
contribution Sc'^^^T), Eq. p.7p . and thus refine this in- 
termediate result. 



Al](2)(T) = -_^ J dndndKdK^/2(q) 

Q 

x{4[7„~j'+[7„-j'tr[a^a']tr[a^a']} 

X {uiyvFin-ci)e + 216*} {-uiyvpin-qje + 2ir} 

X ^^~ui'vp{h-q)e + 2ie*^ |wwi.(n-q)^ + 2ir| 

.gn(g)5fi(-Q) , (3.3) 

where gn{Q) is the bosonic Green's function, Eq. (|2.54p . 
Performing the remaining integrations over the Grass- 
mann variables, evaluating the spin traces, and using 
Eq. p.ip , we obtain the relevant low energy second order 



B. One-loop corrections to 54 

The second order result Eq. p.6p for the thermody- 
namic potential cannot provide the full qualitative pic- 
ture of the non-analytic corrections because logarithmic 
contributions arise in higher orders in the coupling con- 
stants. At sufficiently low temperatures, they become 
large for an arbitrarily weak interaction. In this and the 
next section, we study the logarithmic divergencies in the 
leading one-loop order. 

Considering first the quartic action ^4, Eq. (|2.55p . the 
one-loop order of the expansion in the coupling constants 
yields the diagrams shown in Fig. 01 We want to show 
that only diagram (a) is important and leads to the loga- 
rithmic divergency in the limit T ^ 0, n ^ —h, whereas 
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the contribution of the other diagrams remains finite in 
this limit and does not contain large logarithms. For this 
purpose, we need to focus on the momentum structure 
only while the spin structure present in 54 in the 7*-term 
has nothing to say about the existence of a logarithmic 
divergency. Therefore for the sake of a simpler presenta- 
tion, we only consider the 7'^-term for the moment. 

The logarithmic contributions come from running 
bosonic frequencies with |a;| ^ 2'kT, which are in the fo- 
cus of the following considerations. Contributions from 
the region |a;| < 27rT produce terms of higher order in 7"^ 
and, therefore, only give perturbative corrections to the 
coefficients in front of a large logarithm. 

The standard diagrammatic technique based on the 
contraction rule, Eq. (|2.53p . yields for diagram Fig. |3Ja) 
the one-loop vertex correction 



(a) 



1 f V- 

dndndftdKuu — 

KQ1Q2 



xtr 



xtr 



(3.8) 



written as 



5^c^ JJ!^ [^^^Y 

'nn ^'nnJ / 27r 



./(q_L-qi)/(q2-q_L) 
1 



TV 



VFq\\ + VF{5n-(i±) - ^ci\ 



(3.12) 



As to the radial component g|| . it will be sufficient to 
know the cutoff's order of magnitude while the precise 
form of the cutoff function /(q) is irrelevant to it. As 
a result, it is justified to restrict its dependence to the 
transverse momenta qj^ — as has been done in Eq. (|3.12p 
— while keeping in mind that wf9|| ^ A with A being the 
upper boundary of the bosonic spectrum. 

Integrating over the radial component g|| , we find 



^7-^ = [7'^]2 / ^ 

'nn L /nnJ / 2(Jo 



/(q±-qi)/(q2-q±) (3.13) 



with 



A- 



Auu qo f du! 



uvf tt J27rT 27r uj^ — [vF{Sn-q±)]'^ 



(3.14) 



This vertex SS^""^ reproduces the analytical form of the Remarkably after the integration over q\\ , the 

term q^/(2TO), which is of order A and, thus, in 
principle large, has dropped out. Therefore, the subse- 
quent integration over the frequencies uj leads in the limit 
T, |5n| — to a logarithmic divergency. The important 
observation is that the divergency of diagram Fig. IH^a) 
is not sensitive to the question whether we are in di- 
mension d = 1 or d > 1 and thus the logarithm appears 
in any dimension d. The term VF{Sn-(i±) ~ Wi^gol'^nl, 
though being dependent on the transverse momenta, 
vanishes in the limit i5n — >0 independently from d. This 
however is exactly the limit in which A-^ eventually 

^ nn ^ 

enters physical quantities such as the thermodynamic 
potential Sfl, Eq. ()3.6|) . Explicitly, we find 



quartic action ^4, Eq. (|2.55p . where a correction — (57"— /2 
should now be added to the bare quantity 7'k/(qi— q2). 
This correction — a function of it, u, n, n, qi— q2, and K 
— is determined by the integral over the running four- 
momentum. 
We find 

^7:-, = ^K~,]^E/(q-qO/(q-q2) 
Q 

X 9n+^ [K + Q)g^_2^{K - Q) . (3.9) 

This integral is conveniently calculated introducing the 
angular variables 



n=i(n-n), (3.10) 
Sn = ^(n + n) , 
and their projections of the momentum vector q. 



A — : = 4uu — In — — T- 

Vmax(27rr,VFgo|<5n|) 



with the constant ly* defined as 



1 



qo_ 

2^^ vf 



(3.15) 



(3.16) 



q\\ = (n-q) , 

q± = q 9|n 



(3.11) 



Calculating the integral in Eq. p.9p with logarithmic ac- 
curacy and keeping in mind that the essential K in the 
final integration, e.g. in the diagram in Fig. |3Ib), will be 
of order T, which is the lower cutoff of the logarithms, 
we can safely put = in Eq. p. 91) . Thus, Sj"^ can be 



By construction, the transverse momentum only varies 
on a small arc with a length of order 250 <S pf on the 
Fermi circle, but in final results this arc should extend to 
a semicircle, corresponding to go ~ ij^ I'^)pf or v* ^ v 12. 

Still the integration over q^ in Eq. p.l3p remains to 
be done. Since VFq\\ ~ q5_/(2m) ~ A in the region of 
the logarithm, the transverse momenta are much larger 
than the parallel ones, q_L ^ gy. As this statement re- 
mains true for the leading terms in all logarithmic orders. 
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FIG. 4: One-loop diagrams for the quartic action 54: Diagram (a) is logarithmic independently from the dimensionality of the 
system whereas diagrams (b) and (c) are negligible in dimensions d > 1 due to curvature effects. Diagrams (d) and (e) vanish 
after -integration for any d and, finally, diagram (f) vanishes due to supersymmetry. The contributions of diagrams (a') 
and (a"), each of which is logarithmic on its own, cancel each other. 



we may neglect in the relevant limit of (5n — j-O the paral- 
lel momenta 91,2; || in the cutoff functions in Eq. ()3.13p . 
Then, the remaining integral becomes 



dcy 

2go 



7(q±-qi,-L)/(q2,_L-q±) • 



(3.17) 



This is nothing but a convolution [/*/](q2.±— qi.i), which 
becomes a product after employing a Fourier transforma- 
tion. 



/(r±) = / /(q±) e'^ 



dq_L 
2go 



(3.18) 



The form of the Fourier transform in Eq. p.lSp has been 
chosen such that both r_L and /(rj^) are dimensionless. 
The value rj^ = 1 corresponds to the minimal length of 
the theory, which is given by l/go- In Fourier represen- 
tation, the vertex correction takes the final form 



(3.19) 



with the functio n A^ , Eq. p.lSp . containing the loga- 
rithm. From Eq. p.io)) . we understand that it is actually 
the quantity 7S~,/(rx) which flows during a renormaliza- 
tion procedure. 

In conclusion, the diagram Fig. SKa) logarithmically 
renormalizes the backward scattering amplitude 7^ of 
the quartic action 54 and this logarithmic contribution 
comes independently of the dimension d. Diagram (a) 
corresponds in conventional fermion diagrammatics to a 
rung of the particle-particle ladder. We discuss this cor- 
respondence in Appendix [Cl 

Let us now turn our attention to diagram Fig. |4l[b). 
After the integration over the internal momenta and fre- 
quencies, this diagram also reproduces the structure of 54 
but in contrast to diagram (a), it is not logarithmic for 
dimensions d > 1 . In order to support this statement, we 
consider the vertex correction 8S'^^ , which has the same 
form as 5sf^ , Eq. (jS^, except that A^, Eq. ((XTil) . is 



replaced by the function 



(3.20) 



.9-^ <,! (A" + Q - (Oi + O2)) 



In order to estimate the integral over Q in Eq. p.20p . 
we put all external momenta and frequencies equal to 
zero. Also, the precise form of the cutoff functions /(q) 
is not needed for this estimate and therefore, we do not 
write them here for simplicity. Then, in the frame of 
the angular coordinates from Eqs. p.lOp and p. lip we 
obtain 



(3.21) 



d(7|| dqj 



1 



1 



27r 27r \u:-VFq\\- ^(^^\u:^VFq\\~ ^^^^ 



Similar to the case of diagram Fig. IH^a), the unit vec- 
tors n and n need to be close to anticoUinearity if we 
want to achieve the largest value of the integral. Passing 
from the first to the second line in Eq. (|3.2ip . this has 
already been assumed. Integrating over g|| yields 



A^i^oc 



^E / |,,,| . 1 MB^ 



dqj 



sgnw -2 

1m q± 



27r 



(3.22) 



This intermediate result for A—; already demonstrates 

nn 

what makes diagram Fig. IH^b) essentially different from 
diagram (a): Here, the transverse term q5^/(27Ti) does 
not drop out. Moreover, nothing prevents the momen- 
tum q from being large and the energy q^/(2m) from 
taking values of the order of the cutoff A. As a re- 
sult, logarithms analogous to those that appeared from 
the diagram Fig. |U[a) are suppressed by the presence of 
the transverse term. Since the latter exists for dimen- 
sions c? > 1 only, the vanishing of the divergent contribu- 
tion of diagram (b) is clearly due to the higher dimension- 
ality of the system. In one dimension, diagram (b) would 
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give the same logarithmic contribution as diagram (a). 
Thus, one can interpret the cancellation of the diagram's 
logarithmic contribution as an effect of the finite curva- 
ture of the Fermi surface in d > 1. 

Studying the correspondence of diagram Fig. IH^b) to 
the conventional fermion diagrammatic technique — cf. 
the discussion of Fig. [T^ in Appendix [C] — , it is possible 
to identify this contribution with the particle-hole lad- 
der and polarization bubble diagrams. As is well-known, 
both these contributions are logarithmic in d = 1, but 
not in higher dimensions d > 1. 

Diagrams Fig. HJ^d) and (e) trivially vanish because 
the angular variables n in both the internal propagators 
are necessarily close to each other. Therefore, the in- 
tegration contour for q|| can be closed without residues 
inside and the integral equals zero. Diagram (f) contains 
a closed loop of bosonic propagators. Since this diagram- 
matic substructure is fully supersymmetric in the sense 
of Eq. (|2.48p . the contribution of the entire diagram van- 
ishes. 

For the study of the remaining diagrams, the discus- 
sions of the diagrams Fig. |4|a) and (b) can rather easily 
be extended. The reasons why a logarithmic divergency 
appears in diagrams Fig. Ul^a') and (a") while it is sup- 
pressed in diagram (c) in dimensions d > 1 are similar to 
those used for diagrams (a) and (b). 

Considering the diagrams Fig. |4l[a'), (a"), and (c), one 
encounters an, at first glance, very unpleasant property: 
Neither of them reproduces the formal structure of S4 for 
fixed external momenta. Once more neglecting the cutoff 
hmctions for the moment, we obtain for (a') and (c) the 
analytical expressions 



dndnd/tdii 

KQiQ 



dqi ^bhl' 



X tr 



X tr 



and 



(3.23) 



(c) 



6S, 



X tr 



X tr 



1 

'4^ 



dndndAtdii 



Anqo 2 



KQ1Q2 

, ^(-A' - Qi, (Js: + Q2, k) 



2pp 



(K-Q2, 



2p_F 



(3.24) 



The contributions SS^^"'^ and dS^"' , Eqs. ([3:23)) 
and p.24p . contain the logarithmic function A^, 
Eq. (|3.15|) . This function is sensitive to deviations 
of n from — n with n ~ — n being the only region of im- 



:(c) 



portance in our consideration. The structure of the ac- 

fa') (c) 

tion ^4 is formally not reproduced by and SS^ be- 

cause of the presence of the momentum qj in Eqs. p.23|) 
and (IS211). 

However, one can easily see that the momentum 
enters the vertex corrections ds[^ ^ and SS^^ quite dif- 
fcrcntly. As concerns SS^ , the additional rotation 
of the vectors n and n represented by qj^ does not 
change the direction of n and n with respect to each 
other. In other words, if n = — fi, then n + q^/pF = 
— [n — q^/PF] + C[(qi/PF)'^]- This makes the presence 

fa') 

of the momentum q^ in SS\ unimportant. One should 
simply keep in mind that eventually we are to calculate 
the thermodynamic potential correction S^l. Here, SS^ 
can enter in the leading order correction as in Fig. El^b) 
or be a part of a larger ladder containing other iS4-blocks. 
It will turn out that the terms including q^ drop out af- 
ter the integration over the parallel momenta in the 
additional loops. Formally, it is therefore legitimate in 
such diagrams to simply put q^ = 0, thus making the 
diagram (a') give the same logarithmic contribution as 
(55^ \ Eq. (IMl). 

In contrast to SS[^ ^ , the vertex SS^^ is entered by q^ 
in an asymmetric way changing the mutual direction of n 
and n. As a consequence, if we attach another 54-block to 

(c) 

the right of SS^ , a curvature term of order A containing 
q^ will necessarily survive the parallel momentum inte- 
gration, resulting similarly to the scenario of diagram (b) 
in the cancellation of the otherwise emerging logarithm. 
Inserting SS^'^'^ into the perturbation series for the ther- 
modynamic potential. Fig. [3{b), smears — due to the 
presence of q^ in one of the propagators — the impor- 
tant region around n = — n. The expression of the form 
Eq. (|3.4p in Sec. IIII Al will accordingly be no longer suffi- 
ciently sensitive to the backscattering limit, thus result- 
ing in the suppression of all backscattering logarithms. 
For these reasons, the vertex SS'^ ^ can be excluded from 
the class of the important one-loop diagrams. 

The final one-loop contribution to be dealt with is the 
diagram Fig.[3]Ja"). Its evaluation at small external mo- 
menta yields the exactly same analytical form as SS^^\ 
Eq. (|3.23p . but due to a necessary transposition of the 
Grassmann fields with opposite sign, 



6S, 



(a") 



(a') 



(3.25) 



Thus renormalizing the quartic action using the renor- 
malization group, diagrams (a') and (a") cancel each 
other. One straightforwardly checks that this cancella- 
tion still prevails when we consider both charge and spin 
channel of the quartic action S4, Eq. (|2.55p . at the same 
time. 

As a result, the first-loop analysis of the quartic ac- 
tion 54 clearly demonstrates the existence of logarith- 
mic divergencies arising due to the interaction of the col- 
lective excitations of the Fermi gas. These divergencies 
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(b) 



^ (b)^ 



54 



interaction in Sec. IIII Bl we refrain from an explicit dis- 
cussion. 



FIG. 5: Logarithmic one-loop corrections to the quadratic 
action 52, Eq. (|238|l . 



IV. RENORMALIZATION GROUP 



(a) 



S4 



(a') 



FIG. 6: Diagrams (a) and (a') represent the logarithmic one- 
loop corrections to the cubic action S3, Eq. (|2.57|l . 



originate from just one of the various one-loop diagrams 
shown Fig. 21 namely diagram (a). 



C. One-loop corrections to S2 and 53 

In the previous section, we have identified the relevant 
logarithmic one-loop corrections to the quartic interac- 
tion S4. In a general leading logarithmic diagram of or- 
der ln"(A/T) with n being a large integer, nearly all loga- 
rithmic factors will be due to 4S4iS4-loops. The quadratic 
and cubic parts ^2 and 1S3 of the action, Eqs. (j2.58p 
and (|2.57p . that are needed to break the BRST sym- 
metry in diagrams for a thermodynamic quantity, serve 
as the "abutments" of the big 1S41S4-I00P structure. 

In principle, the analysis of the one-loop diagrams for 
the quadratic and cubic vertices is very similar to that 
performed for 54 in Sec. IIII Bl Let us begin with the 
quadratic action S2- One- loop corrections to 52 come 
from diagrams built from and 52 54. As in the 

renormalization of the quartic action, most diagrams are 
negligible since they vanish due to supersymmetry [like 
diagram Fig. lUJf)], as a result of integration over q|| in 
cases when the integrand is an odd function of g|| [like di- 
agram Fig. Hljd)] , or due to higher-dimensional curvature 
effects as in the case of diagrams Fig. IHb) and (c). The 
only diagrams yielding in fact a logarithmic contribution 
proportional to A^, Eq. p.lSp . in dimensions d > 1 are 
those shown in Fig. [5] 

Corrections to the cubic action 53 appear in form of the 
one- loop diagrams built from 5354. Logarithmic renor- 
malizations come from the diagrams shown in Fig. EJa) 
and (a'). 

One-loop diagrams different from these two do not con- 
tribute logarithmically. This is once more a consequence 
of symmetry aspects and curvature effects. Discussing 
the cubic one-loop vertices, one should also bear in mind 
that they finally affect the thermodynamic potential only 
via the effective quadratic vertex in Fig. [SJa). Since the 
arguments follow the same reasoning as for the quartic 



With the perturbative analysis from the preceding sec- 
tion, all the relevant logarithmic one-loop diagrams are 
at hand and we are ready to apply a one-loop renormal- 
ization group (RG) scheme. At the end of the day, the 
energy scales above T will be integrated out of the field 
theory and we will obtain the specific heat in terms of the 
basic backscattering diagram Fig. E^a) with renormalized 
coupling constants. 

It is important to mention here that the renormaliza- 
tion of the quartic term 54 differs from the renormaliza- 
tions of 53 and 52. The former can be obtained both us- 
ing the RG scheme and summing ladder diagrams, while 
the latter ones do not allow for a study based on simple 
summations of ladder diagrams. 



Generalized action 



The actions 54, 53, and 52, Eqs. (|2.55p ~ (|2.58p . con- 
tain various subterms which in general have a specific 
fiow behavior under the RG action. In order to facili- 
tate the RG procedure, we generalize the action a priori 
and introduce proper coupling constants. As a result, we 
write the quartic interaction as 



54 = — — J dndndftdK uu r'^(q^) 

KKQ 

X tr [*„(-i^, K)*n+JL^ (K + Q, k) 
xtr *fi(-X,K)*^_^(X-0,K) 

/dndnd^dK uuT^(<i±) 



(4.1) 



X tr 
X tr 



KKQ 

*„(-A', n)a^^^^^ {K + Q, k) 

2pp 
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the cubic action in the form 

iSa = ^ / dndnd^dK^M 

^ •' KQ 



(4.2) 



X tr 



+ — [ dndfidftdK ' 



According to Eq. (|4.4p . the cutoff functions /(q±) are 
absorbed into the couphng constants. This is a conve- 
nient definition since we have seen in Eq. p.l9p that the 
quantities flowing with the RG are 7!-C*/(q±) rather than 

the interaction constants 7-^^ themselves. The coupling 

nn 

constants with index "_L" also contain the transverse mo- 
mentum qj^ as prcfactors. As the flowing coupling con- 
stants in Eq. (|4.4p arc functions of qj^ , we arc formally 
applying a functional RG procedure. 



KQ 



X tr 



X \^uiyvF9 (5n-Bi(q^) + 2ir Bg(q^)} tr P*f,(-Q, k)] 
and the quadratic action as 
& = f dndndftdKy^ 



(4.3) 



xtr [*„(Q,K)]tr h«'fi(-Q,K)] 

X |(2ir)(2ir) A^aiqj^) + 2{2ie*){wvF0 Sn) Al(q^) 
+ [uvvpO 6n){uvvp9 5xi) A'|[^(qj^)| 

[ dndnd^dKy^ 

J ^ 



Q 



tr [a'^-^niQ, 1^)] tr [^a'^^^i^Q, ^)] 

{(2ir)(2ir) A-^{^^) + 2{2ie*){ui^VFe Sn) Al(q^) 

+ {uvvpB 5n){uvvp9 Sn) A5^^(qj^)| . 



Equations (j4.ip - (j4.3p reduce to the original formulation, 
Eqs. (|2.55p ~ (|2.58p . if we insert the bare values for the 
coupling constants, which are given by 



r^/^(qj 



Bf (q^)^ = A^^^(q^)^=7:4.V(q±)^ 



c/s 



c/s 



Bf (q^)^=Al'^'(qx) 



c/ s 



A-I(qx) ^ = -ql 7„^7(qx) 
(4.4) 

The notations for the (one-dimensional) angular vari- 
able (5n and the transverse momentum q^ arc taken from 
Eqs. (|3.10p and (|3.11l) . Since the coupling constants even- 
tually enter only at backscattering, nn = tt or (5n — > 0, 
the angular dependence is not written explicitly. Note 
that in Eqs. (j4.ip - (j4.4p . the parallel momenta q\\ in both 
the cutoff functions and the prcfactors have been omit- 
ted. The former is justified according to the discussion 
preceding Eq. p.l7p . In the prcfactors (n-q) of the origi- 
nal actions and ^2, Eqs. (|2.57p and (|2.58p . the parallel 
momenta ^|| are irrelevant for the non-analyticitics fol- 
lowing the discussions after Eq. (jA5l) and in Appendix IB] 
This justifies the latter. 



B. Renormalization group equations 

' We develop an RG scheme using the momentum shell 
integration. In one RG step, the large energy cutoff A 
for the one-dimensional parallel spectrum wf9|| (and the 
Matsubara frequencies) is reduced to a still large but 
much smaller cutoff A' ^ A by integrating out the 
fields with parallel momenta of orders between A'/vp 
and A/vp. This yields an action at the energy scale A' 
with renormalized coupling constants. Repeatedly ap- 
plied RG steps make the coupling constants flow. This 
RG flow stops at the latest as soon as the cutoff ap- 
proaches the order of the temperature T. 

In this work, we study the RG flow of the coupling 
constants in Eq. (j4.4p . which comprise the physics of 
the anomalous low energy behavior of a two-dimensional 
Fermi liquid, at the leading one-loop order. This corre- 
sponds to a summation of all orders of ln(A/T) at leading 

order in 7Jr^*. The relevant one- loop diagrams have been 
completely identified in Sec. IIII Bl 



RG equations for S4, 



Equations (|3.8p . (|3.15p . and p.l9p determine how the 
coupling constant r'^(qj^) in the quartic interaction ^4, 
Eq. (gH), or its Fourier transform r'=(r_L), cf. Eq. ((XT^ . 
is renormalized if the energy cutoff A is reduced to A/5 <C 
A. Defining 



(4.5) 



the correction dr'^(r^) to the coupling constant r'^(rj_) 
in one RG step would be 



1 



dr^(r^) = --[r^(rx)]'de 



(4.6) 



if we could neglect the spin channel r*(rx). 

Including the spin channel, we have to examine the 
spin structure of the one-loop diagram Fig. lUJa). Equa- 
tion (|4.6p is the result of attributing both 54-blocks to F'^. 
Replacing in one of these blocks F'^ by F*, reproduces 
the spin structure of the spin F^-vertex. If both blocks 
belong to the spin channel, the algebra of the Pauli ma- 
trices allocates renormalizations to both the F'^ and the 
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vertices. Then, Eq. (j4.6|) should be replaced by the 
RG equations which are written as 

^^ = -^{[nr.r + 3[r^(r,)]^}, (4.7) 

= -l{2r^(r J r^(r J - 2[r^(r^)]^} . 

This system of two differential equations is decoupled for 
the linear combinations 



ri(ri) =r^(ri)-3r«(ri) 
r"(r^)=r(r^) + P(r^), 



(4.8) 



whose bare values are given by 



r^K) 



./(r±) 
r"(r^) 



37- 



(4.9) 



/(r±) 



7^^ + 7^ ■ 

nn nn 



In terms of r^/^^(rj^) 
form 

drVii(r 



the RG equations (|4.7p take the 



(4.10) 



Equation (|4.10p with Eq. (|4.9p as boundary condition is 
easily solved, yielding 



r"(ri;0 = 



1 + i7^/(r.) ^ 



(4.11) 



where in the relevant backscattering limit 6n 0, the 
quantity ^ varies between and 4:uu{i/* /i>) ln(A/T). 

The renormalized coupling constants T^'^^, Eq. (|4.1ip . 
can also be obtained summing the relevant ladder dia- 
grams built of the quartic vertices. Such a ladder is con- 
structed adding rungs of 54-blocks one by one in the way 
corresponding to the diagram in Fig. |4Ka). Considering 
the correspondence of the boson diagrams to the conven- 
tional fermion ones as discussed in Appendix [Cl one can 
see that these ladders are related to the usual particle- 
particle Cooper ladders. The evaluation of the arising 
geometric series is depicted in form of diagrammatical 
Bcthe-Salpether equations in Fig. [T] The two equations 
decouple in a complete analogy with Eq. (|4.8p and even- 
tually yield the same result Eq. (j4.1ip as obtained using 
the RG equations. 

The correspondence between the ladder form of the 
renormalized quartic vertex and conventional Cooper 
ladders suggests that the logarithmic renormalizations of 
the coupling constants can be attributed to the supercon- 
ducting correlations which in case of an attractive inter- 
action cause a phase transition toward superconductivity 



at a critical temperature Tc. Re-expressing the bare cou- 
pling constants 7^^^^, Eq. (|4.1ip . in terms of the original 
interaction potential, Eq. (j2.56p . we come at backscatter- 
ing nn = TT to the relations 



7; 



7? 



Vi2pF)} 
V{2pf)} 



(4.12) 



In final results, contributions containing can be inter- 
preted in terms of a spin singlet while contributions due 
to 7^^, which enter final results with a prefactor of three, 
can be attributed to spin triplets. Note that the latter 
coupling constant vanishes in models with a contact in- 
teraction. 



2. RG equations for S3 

The relevant one-loop diagrams for the cubic interac- 
tion 53, Eq. (|4.2p . are shown in Fig. IHlJa) and (a'). In 
contrast to the quartic interaction 54, the renormalized 
cubic vertex cannot be obtained from simple ladder sum- 
mations and an RG procedure seems unavoidable. One 
can understand this fact from Fig. [51 which shows that 
one has always two different choices in attaching another 
quartic block when building higher-order logarithmic di- 
agrams. As a result, a diagram for a general leading 
vertex correction acquires a rather complicated topology 
and momentum structure. 

Similarly to the RG equations for ^4, Eq. ()4.7|) . the 
RG equations for the cubic action take the form 



dB;;/^(r x) 



d^ 



= -{r^(r^) BS/^(r^) + 3r-^(r^) Bg/^(r^)} , 

= -{r^(r^) B^/^(r^) + r^(r^) B;5/^(r^) 

-2r^(r^) B;5/^(r^)} . (4.13) 

the linear combina- 



Inserting in analogy with Eq 
tions 



Bi/^(rx) = B;;/^(r^)-3B^/^(r^), 



'0/ 



o/±\ 



(4.14) 



Bii/^(rx) = B;5/^(r^)+Bg/^(r^) 



'0/ 



0/_LV 



into the the RG equations in Eq. (|4.13|) reduces these 
equations to the form 



dB^/" (r , ) 



(4.15) 



With the knowledge of the renormalized coupling con- 
stants ri/"(rx) of the quartic action, Eq. (|ITT|) . the RG 
equation Eq. (|4.15|) is nothing but a homogeneous linear 

differential equation for the coupling constant B^'^^^ (r^ ) . 
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FIG. 7: Diagrammatical Bethe-Salpetlier equations for the <S4-ladders. Dark gray blocks stand for the renormalized vertices 
while light gray blocks stand for the bare quartic vertices. V and are defined to project a general quartic vertex onto its 
vertex component of the form or [*cr''v[>] [*ct''*], respectively. 



Using the boundary conditions 



BY"(r.) 



= 7-4 / r_L) 



(4.16) 



I/II 

= 

5=0 



= 7rs |i9oi9r^/(rj 



cf. Eq. (j4.4p . the solutions of the RG equations in 
Eq. (|4.15p are found to be 



Br(r.;0 = 



7-4 .f(r_L) 



(4.17) 



[1 + z^'-^) ^] 

Equipped with the knowledge about quartic and cubic 
coupling constants at any energy scale, we are now ready 
to finally determine the renormalized coupling constants 
of the quadratic interaction 1S2 . 



3. RG equations for S2 



The quadratic action 1S2, Eq. (j4.3p . is renormalized 
according to the one- loop diagrams in Fig. EJa), (b), 
and (b'). The spin structure of these diagrams is com- 
pletely analogous to the one-loop corrections to the quar- 
tic and cubic terms. Introducing corresponding combi- 
nations of the coupling constants in Eq. (|4.4p as 



^o/_l/_l_l('^J-) — ^0/±/J 
(r^) = 



A" 



^0/±/J 



_(r^)-3Ag/^/^^(r^), 

(4.18) 



we immediately write the decoupled RG equations. 
These are 



dAy"(r^) 
dAr(r,) 



de 

dAYi^(r,) 



de 



ri/"(rj 
ri/"(rj 
ri/"(rj 



AY"(r,)-Br(r.)BY"(r,) 

AYi\r,)+[BY"(r,)]\ 

(4.19) 



Since T^l^\vi_), Eq. (HHI) . and B^/" (r^), Eq. BTf)) . are 
known functions of ^, we are once more to solve linear 
differential equations, which are now — in contrast to the 
RG equations for the cubic action, Eq. (|4.15p — non- 
homogeneous. The boundary conditions for Eq. (|4.19p 
take the form 



A^/"(r,) 
AY"(r,) 
AYi^(r,) 



5=0 
5=0 
5=0 



i/ii 

nn 

I/II 
nn 



7^ /(rj 



= 1^ [igo 9r^/(i-_L)] 



(4.20) 



cf. Eq. (|4.4p . As a result, the solutions of the RG equa- 
tions ()4.19p are 



.1/11/ I/II J./ N 

Aq' (r±;0 =lh. f(^^) 



1 - h^V(r.) ^ 



[1 



AY"(r^; 0=190 



9 



I/II J-/ N 

7-4 /(r± 

nn 



AYi^r,;0 



2 
90 



9rj 
92 



7iX"/(r.) 



9r 



(4.21) 



(4.22) 



(4.23) 



Equations ()4.2ip - ()4.23p constitute the final result for 
the renormalized quadratic vertex and complete the 
RG study of our low energy field theory, Eqs. (j4.ip - 
(|4^ . Stopping the RG flow at ^ = 4uu(i/*/i/) ln(A/T) 
yields the renormalized coupling constants in the relevant 
backscattering limit Sn 0, which enter the effective ac- 
tion after integrating out all superfields at energy scales 
larger than T. Therefore in low temperature calculations 
of thermodynamic quantities based on the renormalized 
action, the summation of all orders in the large loga- 
rithm ln(A/r) is at leading order in j^^^^ automatically 
included. 



V. SPECIFIC HEAT 

As a result of the renormalization group procedure 
developed in the previous section, we have the full 
knowledge about the physics of our low energy theory. 
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Eqs. (|4.ip - (|4.4p . at any energy scale A' < A. The rele- 
vant energy scale for the calculation of the specific heat 
is according to Eqs. ()3.ip and p.2p of the order of the 
temperature T A. So, inserting the quadratic vertex 



with the renormalized coupling constants, Eqs. (|4.2ip - 
(|4.23p . into the formula for the thermodynamic poten- 
tial -(T/2)(5|}, we obtain instead of Eq. dM]) the fol- 
lowing expression: 



SniT)= Ty:-fp] /4^dndn- , ^^^^^^ , 

^ ' y ^ J 2ti j J (27r)2 [iLo-VF{n-q)][-iLO + VF{n-ci)] 

< j'j\uAu uu[{[/^\{cir,i)]^ - /^l{cir,0^\A^^^^ . (5.1) 



For i = Auu{v* /v) ln(A/T), Eq. (|5.1I) yields the anoma- 
lous contribution to the thermodynamic potential in all 
orders in the logarithm ln(A/T) ~ \Ti{£p/T) at leading 

order in the bare couplings yj^^ ■ This means that in 
the limit of low temperatures T <^ ef considered here, 
Eq. (|5.ip gives the full physical result. Explicit results for 
the specific heat can be extracted by expansion or speci- 
fying the cutoff function /(qj.) that enters the renormal- 
ized coupling constants, cf. Eqs. (|4.2ip - (|4.23p . 

An important observation in Eq. (|5.ip is that the en- 
tire contribution is expressed as a sum of two contribu- 
tions. The first one contains only the bare coupling con- 
stant 7^ while the second contribution contains only the 
coupling constant 7^^. In other words, the fluctuations 
attributed to the spin singlet [7^] are in leading order 
completely separated from the fluctuations due to the 
spin triplet [7^^], cf. the discussion in Sees. IV Cl and IV Dl 

A. Low order perturbation theory 

Before evaluating the specific heat in all orders for 
a certain cutoff function /(r^), let us explicitly calcu- 
late the non-analytic contribution to the specific heat 
in the third order in the interaction. This approxi- 
mation is valid for not very low temperatures T, such 
that the quantity ^y^^hi{K/T) is still small, i.e. if 
T > Aexp(-[7y"]-i) is fulfilled. Resuhs in this limit 
are known from conventional perturbation theory2^ and 
therefore, we can check our bosonization approach and 
see how it works. In Appendix |Bl we recalculate the 
third order starting from the bosonic diagrams rather 
than from Eq. (|5.ip . 

In the second order in 7]/^^ the renormalized vertices 
in Eqs. arc reduced to 

Ai'^rJ - 7a"/(r J - 2[iBV(rxf « , (5,2) 

A./"(rJ.i,„^{o^:i"/(r.)-[,lf/(rJ]^«}, 



Returning to the momentum representation by Fourier 
transformation, we find 

A^/"(q^) ^ 7^"./(q±) [./*/](q±) i , (5.3) 

AY"(q^) ^ -q47lX"/(qJ " Wi]' [/*/](q^) , 

with [/*/] denoting a convolution as in Eq. p.l7p . 

Inserting Eq. (j5.3p into Eq. ()5.ip . using the relation ^ = 
Auu{v* /v) ln(A/r), and taking the result from the cal- 
culation in Appendix |^ we find the non-analytic ther- 
modynamic potential correction in the third order in the 
interaction as 

<5f^^^{[7iP + 3[7if 

-^[/*/](0)([7y^ + 3[7if) ln(A)}. (5.4) 

Choosing the cutoff function as /(q±) = 0(go — 
the prefactor [/*/](0) just gives unity. In this case, 
we can also estimate v* /v ~ 1/2, following the lines 
after Eq. p.l6p . Finally applying the thermodynamic 
relation 5c = —Td^6fl/dT^ yields the third order non- 
analytic speciflc heat contribution at low temperatures 

TTVp L 

-2([7if + 3[7y]^) ln(^)}. (5.5) 

This perturbative result is applicable in the limit of not 
very low temperatures such that on one hand ln(A/T) ^ 

1 while on the other hand 7y^^ln(A/r) <C 1. 

Equation (jS.Sp contains the logarithmic contributions 
from both spin and charge excitations, which corrects the 
result of Ref . [l^, where the contribution of the charge ex- 
citations was missed due to the neglect of the curvature 
of the Fermi surface. (However, we emphasize that the 
results of Ref. [l^ remain correct for onc-dimensional sys- 
tems.) At the same time, this correction to the specific 
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heat 5c, Eq. (j5.5|) . fuUy agrees with the result of the later 
publication by Chubukov and MasloAJ^^. 

To be more specific, the result for the third order non- 
analytic contribution to the specific heat at low temper- 
ature by Chubukov and Maslov can be recasted into the 
form 

ZTTVp 

+ 3(uo-iix)(K-w.-^)')} THn{A/T) (5.6) 

Here, is the scattering angle, = uV{2pp sm{i}/2)), 
and (g^) is the angular average for an arbitrary function 
g^. Following the decoupling into soft modes, Eq. (|2.7p . 
angular averages (g^) are to be replaced by (go + <?7r)/2 
in our model. Applying this correspondence to Eq. (|5.6p , 
we recover immediately the same logarithmic dependence 
of the specific heat (Jc'^-* at third order as in Eq. (|5.5p . 

Thus, on one hand Eq. (|5.5p serves as a good check of 
our low energy model, on the other hand we confirm the 
estimate i'* ^ vjl discussed after Eq. (|3.16p . 

In the remaining of this analysis, we extend the pertur- 
bative result (j5.5p by including the leading in ^-1'^'^ terms 
of all orders in the logarithm ln(A/T). This calculation 
shall complete the picture of the non-analytic corrections 
to the specific heat at low temperatures T. 



B. Full low temperature result 

In this section, we will extract from Eq. (|5.ip the 
anomalous contribution to the specific heat in all orders 
in ln(A/T). As a result, we obtain the full picture of the 
non-analyticities in the Fermi liquid thermodynamics at 
low temperatures T. 

In order to accomplish this task, we should choose a 
model for the cutoff function /(q±), which controls the 
two soft modes represented by Figs. [TJa) and (b). Fol- 
lowing Ref. [isl a suitable candidate is a Lorentzian of 
the form 



/(qx) = 



'-1^1 



1 ■ 



(5.7) 



This choice implies /(O) = 1 so that the result from sec- 
ond order perturbation theory, Eq. p.6|) . remains the 
same. We recall that within the low energy theory, it is 
implied that go ^ Pf- Fourier transforming Eq. (|5.7p 
yields 



Eq. (|5.ip . Since only the leading quadratic in term 
of the expression [aY"]2 - Al/'^Afl is relevant, we may 
neglect higher orders in qj^ from the beginning. Thus, 
the evaluation of the Fourier integrals yields 



aI/II/- n 



i/ii 

nn 



,1/11 



(q±) 



-q± 



I/II ^-^2 ' 
nn 

I/II 

nn 

1 + f 7^X" e 

^ nn 



^ ^ nn ^ 

7. 



4 



A/x(q±) = -ql — ln(l + -7J, ^ 



(5.9) 

(5.10) 
(5.11) 



Applying ^ = 4uu(:^*/zy) ln(A/T) to the renormal- 
ized couplings, Eqs. (|5.9p - (|5.1ip . inserting them into 
Eq. (|5.ip . using the integral 



1 pi 



dudu 



uu + x ^ \n{l -\- xuu) ln^(l-|-a;) 



Jo 



(1 + xuiif 



2x2 



and adopting the result of Appendix [X] for the remain- 
ing integrations, we obtain for the low temperature non- 
analytic part of the thermodynamic potential 5^1 the for- 
mula 



5n = C(3)T3 r ln^(l+^ ^ ^ ij,2(^ ^ ^11^) ■ 



TTUr 



L2 



L2 



Herein, the quantity L is defined as 

'A 



L = ^ In^ 



(5.12) 



(5.13) 



with u* given by Eq. p.l6p . The bare coupling con- 
stants 7y^^ are expressed in terms of the original fermion 
interaction potential V as 

7y" = ^{nO)±n2pF)} , (5.14) 



cf. Eq. KV2\ . 

Equation (|5.12p constitutes our final result for the non- 
analyticities of a two-dimensional Fermi gas with repul- 
sive interaction. In the following, we discuss the correc- 
tions to the specific heat of the Fermi liquid and possible 
instabilities. 



/(rj 



(5.8) 



according to the definition of the Fourier transform in 
Eq. (j3l8| . 

Inserting the model function /(rj_), Eq. (|5.8p . into 
the renormalized quadratic vertices, Eqs. (|4.2ip - (|4.23p . 
we are in a position to Fourier transform them to the 
momentum representation, which is needed for formula 



C. Corrections to the Fermi liquid 

For the Fermi liquid model, the thermodynamic poten- 
tial correction (SO, Eq. (|5.12p . is a regular function for all 

relevant 7^^^^^. By means of the formula 



(5c = -T 



(5.15) 
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FIG. 8: The non-analytic correction to the specific heat of a 
two-dimensional Fermi liquid for a repulsive contact interac- 
tion as a function of •yl^L. The logarithmic-linear plot in the 
background illustrates the behavior at large 7^1/, i.e. in the 
asymptotic limit of very low temperatures T. 



we find the anomalous contribution to the specific heat 
at low temperatures T in the form 



Sc- 



nvr 



£2 



£2 



(5.16) 



This formula has been obtained from Eqs. (j5.12|) 
and (|5.15p by differentiating only the T"^ prefactor but 
not the quantity L, Eq. (|5.13p . Derivatives of L would 
yield subleading terms in the most interesting low tem- 
perature limit. 

Equation ()5.16p constitutes our main result for the 
anomalous correction Sc to the specific heat of the two- 
dimensional Fermi liquid at low temperature in the weak 
coupling regime. Figure [8] shows the plot of dc divided 
by the second order result p.7p as a function of 7^L. 
For simplicity, only the case 7^^ = (contact interaction) 
is plotted. For a general a contribution of the same 
shape as in Fig.[5]but dependent on j^L has to be added. 

With a further decrease of the temperature, the quan- 
tity L, Eq. (|5.13|) . grows such that the non-analytic 
part Sc of the specific heat of the Fermi liquid asymp- 
totically reaches zero as — (In^L)/^^. This asymptotic 
behavior fully agrees with the estimate by Chubukov and 
Maslo\^ in the limit of small z/*/i/. This estimate was 
based on the conjectured relatio n^2,i3 



Sc 



3C(3)T2 



[/f(7r) + 3/,2(7r)] 



(5.17) 



where /c(7r) and /s(7r) denote the charge and spin 
components of the fully renormalized backscattering 
amplitude.— The asymptotic agreement between our 
bosonization approach and the conjecture ()5.17p sup- 
ports its validity. Moreover, full agreement between our 
result Eq. (j5.16p and Eq. (|5.17|) can be shown in the 



limit V* jv <^ 1 in all orders in L for a long-range inter- 
action, considering the potential to be non-zero only 
close to q = 0.— We also note that in the supersymmet- 
ric approach of Ref. [l^ where certain effects of the Fermi 
surface curvature were neglected, the dependence on L of 
the spin contribution to the anomalous specific heat was 
found to have the same analytic shape as the result of 
the present work for the contribution due to both charge 
and spin excitations, Eq. (j5.16p . 

At not very low temperatures, 7y^^ln(A/T) ^ 1, the 
full result for the anomalous specific heat, Eq. (|5.16p . re- 
duces to the perturbative result of Eq. (|5.5|) [except for 
a prefactor of 7r/4 in front of the third order term which 
is specific to the choice of the cutoff function, Eq. (|5.7p ] . 
The logarithmic terms in Eq. (|5.5p were attributed to 
the contribution of superconducting singlet and triplet 
fluctuations (fluctuation contribution of s- and p-wave 
Cooper pairs) <^ This is supported by the form of the 
coupling constants 7^ and 7^^, see Eq. (|5.14p . The com- 
bination of the interaction amplitudes V^(0) and V(2-pp) 
precisely corresponds to what one obtains when summing 
up the singlet [for 7]^] and triplet [for 7^^] Cooper ladders. 

Following this line of reasoning, we interpret the first 
and the second terms in Eq. (j5.16p as contributions 
coming from singlet and triplet superconducting fluctu- 
ations. Of course, for models with a contact interaction, 
V = const, the coupling constant 7^ vanishes and only 
the singlet Cooper pairs contribute. 

It appears that flnding the anomalous correction to 
the specific heat 5c, Eq. (|5.16p . from conventional dia- 
grammatic expansions is rather difficult. Identifying the 
Cooper "wheels" as the relevant diagrams and especially 
identifying the soft modes in these diagrams for an ar- 
bitrary order in the perturbation theory can be rather 
tricky, which is why earlier works^^ restricted themselves 
to low orders while estimating infinite order results by 
plausibility. The bosonization approach presented in this 
paper allows for a field theory study based on simple el- 
ementary diagrams that can be used as building blocks 
for a subsequent rcnormalization group analysis. As a 
result, an explicit derivation of the low temperature non- 
analytic specific heat of higher-dimensional Fermi liq- 
uids becomes possible in all orders in the large loga- 
rithm ln(A/T). 

Formula (|5.16p for the anomalous specific heat con- 
tribution Sc has to replace the well-known second order 
contribution Sc^"^^ ^ Eq. (|3.7p . as soon as the logarithm 

in L, Eq. (|5.13p . becomes of the order of [7^^^]"^. The 
experimental study on '^Hc fluid monolayers presented in 
Ref. conflrmed the validity of the second order per- 
turbation theory down to temperatures of order 1 mK. 
In this experimental setting, however, the coupling con- 
stant 77r is of order unity, which is beyond the applica- 
bility of our theory as it is based on the assumption of 
weak interaction. Nevertheless, provided the logarith- 
mic rcnormalization remains valid at least qualitatively 
also for strong interactions, the logarithms should be de- 
tectable in measurements, although one might need to 
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investigate a broader temperature interval including tem- 
peratures considerably below the mK scale. 



D. Instabilities 



According to the main result, Eq. (j5.16p . discussed in 
the last section, the non-analytic corrections to the Fermi 
liquid thermodynamics at low temperatures are small 
as the function 6c{T)/T'^ decays logarithmically in the 
limit of T — >■ 0. This statement remains valid as long as 
both singlet and triplet constants, 7^ cx ^(0) -I- V{2pf) 
and 7^^ oc V^(0) — V{2pf), are positive. The assumption 
of their positivity guarantees that the system is in the 
Fermi liquid regime. 

If one of these constants is negative, the function ln(l-|- 
^y^^L) in the thermodynamic potential Eq. (|5.12p . 
approaches a pole at |7y^^|_L = 1 when lowering the tem- 
perature T, i.e. when boosting the quantity L, Eq. (|5.13p . 
This pole corresponds to the divergency of the geomet- 
ric series in the summation of the ladders in Fig. [7] for 
^y^^L = — latw = M = l and r_L = 0. The emergence of 
this singularity means the breakdown of the perturbative 
approach close to and beyond |7y^^|L = 1. 

For a repulsive interaction, the coupling constant 7^ 
always remains positive. It would be negative if the ef- 
fective interaction V were negative, corresponding to at- 
tractive interaction. We can identify the resulting sin- 
gularity in the thermodynamic potential with the well- 
known Cooper instability. As a result, one obtains the 
conventional s-wave superconductivity below a certain 
transition temperature Tc obtained from the condition 

So, while 7^ is strictly positive for a repulsive interac- 
tion, the triplet coupling constant 7^^ becomes negative 
as soon as V{2pf) > V{0)- This may actually happen 
for a repulsive interaction. For instance, for contact in- 
teraction, V{2pf) = V{Q) and already a small increase of 
V{2pf) — for instance due to the closeness to a quantum 
critical point — may make the constant j]^ negative. If 
this happens, one comes again to an instability in the 
thermodynamic potential Sfl, Eq. (|5.12p . at j^L ~ —1. 
This instability should correspond to the triplet p-wave 
superconducting pairing. One can come to this conclu- 
sion recalling once more the superconducting particle- 
particle ladders in the spin triplet representation. It is 
easy to see that for a small angle and backward scattering 
the combination V{0) — V{2pf) enters the pre- factor in 
front of the logarithms, which confirms our assumption. 

The condition for criticality j^^L = — 1 determines the 
critical temperature Tc, 



A 



exp 



( ^\ 



(5.18) 



corresponding to the transition temperature into the 
triplet superconducting state. The nature of this transi- 
tion is similar to the Kohn-Luttinger transition towards 



p-wave pairing^ although in the Kohn-Luttinger sce- 
nario, the inverse coupling constant l/|7]r^| enters the ex- 
ponent of the critical temperature Tc, Eq. (|5.18|) . with a 
different power 1^1^ 

Everywhere in this paper, we considered the limit of 
weak interactions. A more interesting situation may oc- 
cur near a quantum critical point (QCP) of a transition 
into a magnetic or charge density wave state. In the 
vicinity of such a point, the collective mode propagator 
X (q, ^) dressed by the electron-hole bubble can be writ- 
ten as 



x(q,f^) 



Xo 



7(|f^l/|q|)' 



(5.19) 



where xo is the susceptibility, determines the close- 
ness to QCP [^~^ = at QCP], and 7 is proportional 
to the interaction. The propagator x(q, i^), Eq- (|5.19p . 
is written in the limit of small bosonic Matsubara fre- 
quencies Q. Higher-order corrections to the propaga- 
tor X (qi ^) can be considered but these are highly non- 
trivial [see, e.g. Refs. [35l - [39| . which invalidates the early 
conjecturc4Si^ that one can describe a QCP by a conven- 
tional 0*-ficld theory with the bare propagator x(q, f^), 

Eq. auim . 

Using the conventional fermionic diagrammatic tech- 
nique for studying the critical behavior near a QCP is 
very difficult. At the same time, the bare bosonic prop- 
agator gn{K), Eq. (|2.53p . in our bosonization approach 
describes directly the electron-hole excitations. Near a 
QCP, it could be modified and take a form similar to the 
one of Eq. (|5.19p . Then, we would be able to derive a su- 
perfield theory with a modified bare action. Within such 
a theory, diagrams to be disregarded for the Fermi liquid 
like, e.g., diagram Fig.|4l[b) are not necessarily small and 
require a special care. This adds to the complexity of the 
theory with very intriguing consequences. We hope that 
our bosonization method can help in studying the QCP 
problem. 



VI. CONCLUSION 

Singling out the low energy spin and charge excitations 
of a higher-dimensional clean Fermi gas with a repulsive 
interaction, we have modified our general bosonization 



scher 



^23.24 



to derive an effective low energy superfield 



theory. This representation allows to conveniently cal- 
culate thermodynamic quantities in the low temperature 
limit. During the derivation, special care has been given 
to the role of Fermi surface geometry in higher dimen- 
sions. As a result, all curvature effects are preserved in 
the final action, correcting the earlier supersymmetric 
approach from Ref. |T^. 

The superfields in our low energy theory are anticom- 
muting but periodic in imaginary time. Consequently, 
the described excitations obey Bose statistics. These 
bosonic excitations include both spin and charge exci- 
tations, interacting in a non-trivial way as described by 
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quartic, cubic, and quadratic terms in the action of the 
superfield theory. 

A perturbative study of the low energy superfield the- 
ory in the backscattering limit yields the well-known 
leading non-analyticities in the thermodynamics and also 
the logarithmic corrections in any dimension d. In di- 
mensions d > 1 however, the class of diagrams producing 
logarithms is narrower than in one dimension as well as 
in the quasi-one-dimensional approximation of Ref. fill , 
where effects of the curvature of the Fermi surface were 
neglected. As a result, the renormalization of the various 
coupling constants in the theory is for d > 1 significantly 
different from the one-dimensional scenario. This is re- 
flected by different renormalization group equations and 
the observation that in higher dimensions, the renormal- 
ization of the quartic part of the action can be understood 
also in the framework of ladder diagram summations. 

The application of the low energy bosonization ap- 
proach to the two-dimensional Fermi liquid and the sub- 
sequent renormalization group analysis have yielded an 
explicit formula for the non-analytic contribution to the 
specific heat 5c, Eq. (j5.16p . This result is of infinite order 
in the large logarithm \-d.{£f/T) and leading order in the 
coupling constants of the weak interaction. As such, it 
is valid for an arbitrary low temperature T ^ ep- The 
dependence of Sc on the logarithm In^ep/T) - plotted in 
Fig. [5]- indicates that the function 5c{T)/T'^ decays as 
1/ ln^(eF/T) for T ^ 0. As discussed in Sec.|VCl our re- 
sult is in full agreement with asymptotic results of earlier 
works based on conventional diagrammatic expansions. 

Remarkably, the thermodynamic potential and the 
specific heat correction, Eqs. (|5.12p and (|5.16p . consist of 
two separate terms of an identical analytical form con- 
trolled by two different coupling constants 7^ and 7]^^. 
The contribution of the term containing only the cou- 
pling constant 7^^ comes with a factor of three as com- 
pared to the term with the coupling 7^. In the discussion 
in Sec. IV CI we interpret these two contributions as com- 
ing from superconducting fluctuations of spin singlet and 
spin triplet types. This statement is supported by com- 
parison of the conventional Cooper ladders with the lad- 
ders of the quartic bosonic action and by the analytical 
form of the effective coupling constants 7^^, Eq. (|5.14p . 
For a contact interaction, 7^^ = and only the singlet 
superconducting fluctuations contribute. We note that 
within our approximation of small fluctuations at back- 
ward scattering, we cannot distinguish between angular 
harmonics of the same parity, e.g. between s- and d- 
pairing, merely distinguishing between the singlet and 
triplet excitations. 

If one of the coupling constants — 7^ or 7^^ — becomes 
negative, the Fermi liquid picture breaks down at a crit- 
ical temperature Tc and a superconducting phase transi- 
tion takes place. The constant 7^ can be negative only 
for an attractive interaction, leading to the conventional 
spin singlet superconducting transition. In contrast, 7^^ 
can become negative also for certain models of repulsive 
interaction. This scenario of triplet superconductivity. 



which is similar to the Kohn-Luttinger one, is discussed 
in Sec. IV Dl Since our perturbative approach breaks 
down close above the critical temperature Tc, Eq. (|5.I8|) . 
the critical behavior itself should be studied introducing 
the superconducting order parameter. 

A very interesting situation may arise near a quantum 
phase transition into, e.g., a ferromagnetic state. Near 
this point, both superconducting and paramagnetic ex- 
citations are important.—"— In the language of the su- 
perfield theory developed here, the non-interacting para- 
magnetic excitations should be described by the bare ac- 
tion, Eq. (|2.51|) . while the superconducting fluctuations 
appear as a result of the interaction between these ex- 
citations. Due to the special form of these excitations, 
cf. Eq. (|5.I9p . the perturbation theory is more compli- 
cated than the one considered here for the Fermi liquid 
and logarithmic contributions may arise in more classes 
of one-loop diagrams. Since calculations in this inter- 
esting situation are not simple within the conventional 
diagrammatic approaches'^"—, we hope that the present 
bosonization technique will become a helpful analytical 
tool for future studies on this topic. 
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Appendix A: Evaluation of the integral in Eq. (|3.4p 

In this appendix, we explicitly evaluate the second or- 
der contribution 6Vt^'^\T), Eq. (lOj) . 

[\uj — vp{n-(i)\[—iuj + VF\n.-<\)\ 

The explicit form of the cutoff function /(q), introduced 
in Eq. (|2.8p . is not important for the second order loop. 
It is sufficient to know that the estimate |q| ^ go *C pf 
holds. Eventually, the calculation shows that effectively 
only those q enter /(q) which satisfy |q| ^ q^. Conse- 
quently, we can effectively put q = in /(q). 

First, we perform the sum and integral over the 
frequency u. The coth- functions resulting from the 
finite temperature Matsubara summation are conve- 
niently combined with the sign functions from the zero- 
temperature frequency integration using the relation 
cothx = sgnx[l -I- 2 ^^-^ exp(— 2/|a;|)]. As a result, we 
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obtain 



WF y (27r)2 'nn (n-q) - (n-q) 



X! {sgn(n-q) cxp 



lvF\{n-q)\ 



T 



sgn(n-q) cxp 



lvp\{h-ci)\ 
T 



)}• (A2) 



From Eq. (|A2p . we understand that relevant "parallel" 
momenta (n-q) or (ii-q) are of order T/vp- It is this 
observation that makes our entire low energy approach 
useful. 

In order to carry out the remaining integrations, wc 
should transform the variables into a frame that would 
better reflect the physics of the scattering processes under 
consideration. For this purpose, we introduce the angular 
variables 

n=^(n-fi), 5n=i(n + n) (A3) 

and the corresponding projections of the momentum vec- 
tor q, 



9|| = (n-q) 



(A4) 



The choice of the angular coordinates and the notation in 
Eq. (jA3[) is motivated by the observation that the most 
important contributions to (5ri'^2^(T), Eq. (jAip . will come 
from small Sn. The phase space region around Sn = 
corresponds to backscattering, nn ^ it, and is also ex- 
actly the region contributing to the logarithmic renor- 
malizations of the scattering amplitude 7^^. 

Transforming the variables according to Eqs. (jA3P 
and (jA4p has to be accompanied by the proper change in 
the integration measure as dndn = 2d(5ndn. Simplifying 
Eq. (|X2|) with the help of the variables from Eqs. (|X3)) 
and (IA4I). we obtain 



d(5ndn— ^ 7?^ 

' nn Z ^ 



VF 



2|9||l 



l9|| 



(<5n-qj.)2 



9 



27r 27r 



exp 



1=1 

T 



(A5) 



The first term 2|g||| in the curly brackets produces an 
analytic and insensitive to small |(5n| contribution to the 
thermodynamic potential. Focussing on the non-analytic 
contributions, we can neglect such terms. One can ob- 
serve at this point that the non-analytic contributions 
arise completely from the term ((Jn-qx)^ in the numer- 
ator (n-q)(n-q) = (^n-qj,)^ — of Eq. (jAip . In other 
words, the momentum components q\\ of the terms with 
(n-q)6' and (n-q)0 in the actions ^2 and ^3, Eqs. (|2.58p 
and (|2.57p . are effectively irrelevant for the analysis of 
the non-analyticities. 

By shifting the integration contour of the parallel 
momentum q\\ as q\\ 1— > iK||, Eq. (jA5[) is reduced to 



Eq. (7.17b) of Ref. 111. Explicitly, Eq. 
casted into the form 



5]) can be re- 



5ili^){T) = - — 

VF 



d(5ndn 



dq± 
27r 



00 

nn Z ^ 

1=1 



X Re 



dKii 



27r k\ + (5n-q^)- 



■ exp 



iIvfh\\ 
T 



(A6) 



The integration over n is trivial, yielding 1 according 
to Eq. (|2.28p . For the integration over i5n, we no- 
tice that the most important contributions come from 
|((5n-qj^)|/K|| < 1 and, since K|| ~ T/vf and q_L ^ 90, 
therefore from the region of small 5w where the esti- 
mate \6r\.\ < T/{vFqo) ^ 1 is valid. This also implies 
that possible logarithms in T of the scattering ampli- 
tude 7^, which we discuss in Sees. IIIIBI and IIII CI be- 
come active. Before we come to the (5n-integral, we inte- 
grate over the "imaginary" momentum Ky and in order 
to facilitate that, we recast the rational integrand as a 
Fourier integral. 



(<5n.qi)2 



1 



2qo\Sn\K 



\/{qo\Sr. 



g-ir_Lqi/90dr_L . (A7) 



Because of the cutoff qo for the momentum qj^, typi- 
cal |r_L| are of order of unity. Inserting Eq. (|A7p into the 
expression for (5f2(2^(r), Eq. (|A6p . we see that the pre- 
exponential is just a power of K|| and the corresponding 
integration is easily performed. 



X ^Re 

1=1 



27^r3 



-ir±ci±/qo 



/ VFqo\Sn\ 
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-3 



VFqod5n dqi 
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According to the estimates discussed in the pre- 
ceding text, essential values of the quantity = 
tiF'?o|'5n|/(r|r_L|) are of order 1. Transforming from the 
integration variable 6n to 0, the upper limit of integra- 
tion will be of order vpqo/T 3> 1, allowing to extend the 
domain of integration to oo. Using the integral 



Re 



(1 



2 



we obtain 
^17(^)(T).^^ 



1=1 



dq_L 
" 27r(jo 



(A9) 



The sum over / just gives Apery's constant C(3), the re- 
maining integrals are trivial, and we obtain for Sn'-^^ (T) 
the result presented in Eq. p.6|) . 
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(a) 
(b) 
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FIG. 9: Diagrams for the third order backscattering contri- 
bution to the thermodynamic potential. 



Appendix B: Third order correction from bosonic 
diagrams 

In this appendix, wc calculate in the leading logarith- 
mic order the anomalous specific heat in third order in 
the interaction by the explicit evaluation of the bosonic 
diagrams. 

As a result of the analysis in Sec. IIIIl the relevant di- 
agrams are those shown in Fig. [Sj Other diagrams arc 
either insensitive to the backscattering region or their 
logarithmic divergency is prohibited by the effects of the 
curvature of the Fermi surface. Since for the diagrams in 
Figs, (nja) and (b) the curvature terms are not relevant, 
we will omit them in the following formulas. 

In analytical terms, diagram Fig. |9l[a) yields the con- 
tribution 



Ar!(3a)(r) = -— f dndn^ 

QQ' 



x([7^]' + 3[7i!.]3) /(q)/(q')/(q + q') 
X {v|(n-q)(n-q) - u|(n-q')(n-q)} 

X 9n{Q)gn{~Q)9n{Q')gn{-Q') (Bl) 

and the second diagram Fig. [SJb) corresponds to 
Af](3b)(r) = _|. J dndii^ 



x([7^]' + 3[7^j3) /(q)/(q')/(q + q') 

X {-4(n-q)(n-q) - w|,(n-q)(n-q') + 24(n-q')(fi-q')} 

X 9n{Q + Q')gn{-Q)gn{Q')gn{-Q') ■ (B2) 

The coupling constants 7!— and j^-K have been intro- 

nn nn 

duced in Eq. (|4J|) . 

The third order diagrams for Arj'^a) (j^) ^nd AfJ^^b) (y) 
contain two loops with correspondingly two running four- 
momenta Q and Q' . Following the idea of Eq. wc 
should subtract the contribution at T = and deal with 
the quantity Sfl'-^^T) = Afl^^^T) ~ An'-^^O) rather than 
with Af7(3)(T) itself. Therefore, one four-momentum ef- 
fectively varies on the scale < T while the other one 
conversely needs to vary on large scales ^ T in order to 
produce the leading logarithmic correction. 



Let us begin the explicit evaluation with the expres- 
sion for diagram Fig.[5l^a), Eq. (|Bip . The first term in the 
curly brackets behaves in a considerably different way for 
the two cases of small or large Q — corresponding to large 
or small Q', respectively. In the case of small Q, the inte- 
gral over Q is calculated analogously to the second order 
integral Eq. p. 41) while the integral over Q' is essentially 
the logarithmic one- loop integral from Ecj^. p.9p . As a re- 
sult, we obtain a correction of relative order 7^^^ ln(A/T) 
to the second order result, Eq. (|3.6|) . 

The opposite case of large Q constitutes an unpleasant 
divergency at large vfQw ■, which is only formally cut by 
the cutoff functions /(q). Fortunately, this ultraviolet di- 
vergency is exactly compensated by an ultraviolet diver- 
gency appearing with the opposite sign in the first term in 
the curly brackets in the expression for diagram Fig.[3t^b), 
Eq. (jB2[) . such that the result is eventually regular. We 
note once more in this context that the parallel momen- 
tum g|[ of the (n • q)-terms in the actions 1S3 and 1S2 is 
irrelevant for thermodynamic quantities, cf. the discus- 
sion after Eq. (|A5|) . Furthermore, since the seeming ul- 
traviolet divergency in vpqw is compensated — an obser- 
vation that is easily generalized to diagrams of arbitrary 
order — , it is completely safe to neglect the part of 
the (n ■ q)-terms in the cubic and quadratic parts of the 
interaction as done in Eqs. (|4.ip - (|4.3p . 

The second term in the curly brackets of Eq. (jBip is 
odd in both Q and Q' and for this reason, one might be 
tempted to disregard that term. However, because of the 
presence of the cutoff functions, the overall integrand is 
not odd in the perpendicular momenta qj^ and q'^, cf. 
Eq. (|3.11l) for the notation. Explicitly, since 



ql /(ql)/(q±+q' 



dql 



q_L 
2 



2go 

/(ql)/(qi 



dq^ 

2(?o 



(B3) 



we observe that the second term of Eq. (|B1[) gives the 
same contribution as the regular part of the first term — 
with one half coming from small Q and one half coming 
from small Q' . 

Now, let us turn our attention to the expression 
for Af](3b)^ (jB2]). The presence of Q' in three de- 
nominators implies that Q is necessarily the large four- 
momentum. The first term in the curly brackets of 
Eq. (jB2[) consequently does nothing more than neutralize 
the ultraviolet divergency in AVl^'^^^ as discussed above. 
The second term is treated in complete analogy with the 
second term in Eq. (|Bip while, finally, the third term 
is effectively of the same form as the first term of the 
diagram Fig.inja) in the limit of small Q. 

As to the cutoff functions /(q), they have played an 
important role in understanding the seemingly odd terms 
in Eqs. (|BT|) and ((B2|) . After applying Eq. ((B3| . all rele- 
vant terms have the form of the first term in Eq. (|Bip at 
small Q. The remaining integration of Q is completely 
equivalent to the second order integral presented in Ap- 



TABLE I: Vertices of the bosonic excitations and related di- 
agrammatical structures of the Hartree and Fock soft modes 
in conventional fcrmion language. 



bosonic 


conventional 
Hartree Focl< 
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pendixl^l There, we learned that the transverse momen- 
tum of the small four-momentum could be safely put to 
zero in the cutoff functions. Thus, the integral of the 
"large" transverse momentum over the cutoff functions 
yields a prefactor of /(O) / /^(q_L)[dq^/2(7o]. For the 
choice /(qj.) = ©(qo — IQ-lI), this prefactor is just unity. 

Collecting all the terms, we find that diagram Fig.[9lja) 
gives 4/9 and diagram (b) 5/9 of the third order correc- 
tion (5ri'^'^'' (r) to the thermodynamic potential, which can 
be written as 

This result clearly agrees with the one obtained from the 
low order expansion of the renormalized coupling con- 
stants, Eq. (|5.4p . 

Appendix C: Boson model versus fermion picture 

In Sec. Ill of Ref. it was checked in the second order 
in the interaction that our method of bosonization allows 
for, in principle, an exact reformulation of the original 
fermion model in terms of bosonic excitations and repro- 
duces exactly each single contribution from the fermionic 
diagrammatics. The choice of a proper diagrammatical 
representation allowed us to identify the bosonic contri- 
butions with the ones of the fermionic picture already on 
the level of diagrams — before explicitly evaluating the 
analytical expressions. 

Following the decoupling into soft modes in Sec. HIl 
the vertices in the bosonic theory collect at the same 
time the fermionic Hartree vertices with a momentum 
transfer close to zero and Fock vertices transferring mo- 
menta of order 2pp. As discussed in Ref. [U, where an 
exact Hartree-like decoupling scheme has been applied, 
the bosonic propagator corresponds to the propagation of 



27 




FIG. 10: Correspondence between the bosonic low energy 
representation and conventional diagrammatics: (a) the bare 
backscattering diagram from Fig. [SJa), (b) its representation 
in Hartree and Fock soft modes according to Table U and 
finally (c) the conventional fermionic diagrams. 

a particle-hole pair in the fermion picture, which will be 
reflected diagrammatically by opposite oriented double- 
lines. Following the derivation of the exact supersymmet- 
ric representation in Ref. [U for eachwise the Hartree and 
the Fock decoupling schemes, we obtain the diagrammat- 
ical representation in Table HI Table U can be understood 
as a dictionary translating diagrams in the boson picture 
into corresponding standard fermionic diagrams. 

As an example, let us consider the diagram for the 
bare anomalous contribution. Fig. fTUT aV Its evaluation 
in Sec. IIII Al returns Eq. p.6p for the correction to fl, 
which yields the leading anomalous T^-term in the spe- 
cific heat c. 

Redrawing the bosonic quadratic vertices with the help 
of Table |T] in all possible ways that the soft Hartree and 
Fock vertices may enter, we obtain the diagrams shown 
in Fig.fTUfb'). Finally, "literally" interpreting the bosonic 
propagators as pairs of opposite directed fcrmion ones, 
we identify the corresponding conventional diagrams, 
Fig.fTUTc). which share the same low energy physical con- 
tent with Fig.lTUTa). Indeed, standard fermion perturba- 
tion theoryi^ yields exactly Eq. p.6p as the anomalous 
contribution to fi, which exclusively comes from the con- 
ventional second order diagrams in Fig. llOf c) . 

One-loop diagrams in the fermion picture 

Figures [TT] and [T^] constitute the correspondence be- 
tween the one-loop quartic vertex corrections 65^^"^ 
and SS^^''' , Fig. |H[a) and (b) , and their contribution in 
the fermionic picture according to Table HI 

Two out of the four fermion lines are seemingly free 
while the remaining two lines interact with an effective 
renormalized interaction. Figure [TT] shows that SSj^^^ cor- 
responds to a particle-particle ladder in conventional dia- 
grams. Particle-particle ladders are known to give rise to 
logarithmic divergencies independently from the dimen- 
sion d, which is in agreement with the analytical result 
for SS^l\ Eqs. (03), (15TB . and (|5TB . As to Fig. [H 
depending on whether at least one of the two iS4-blocks 
is recasted into a Hartree vertex or not, this renormaliza- 
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tion appears in form of either a particle-hole ladder or of 
a polarization bubble. Both of them share the feature of 
being logarithmically divergent in d = 1 dimension but 
convergent in d > 1. Thus once more, the graphical cor- 
respondence reflects the analytical result we have found 
in the discussion of the bosonic diagram, namely the sup- 
pression of the logarithmic divergency due to curvature 
effects in higher dimensions, cf. Eq. p.22p . 



FIG. 11: The bosonic one- loop diagram Fig. [4ja) and its re- 
lation to conventional diagrams. The latter form particle- 
particle ladders, giving logarithms independently from the 
dimension d of the system. 
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FIG. 12: The bosonic one-loop diagram Fig. IHb) and its re- 
lated conventional diagrams. The latter are of particle-hole 
ladder and polarization bubble type, reflecting the absence of 
a logarithmic divergency in d > 1 dimensions. 
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